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ABSTRACT 

The  dynamic  stability  analysis  of  the  airframe  for  consideration 
in  the  control  problem  was  made  considering  six  degrees  of  freedom. 
The  intention  was  to  explore  a  method  and  determine  the  extent  that 
such  a  scheme  could  be  used  in  a  system  known  to  be  non-linear  at 
larger  disturbances.    It  was  determined  that  a  very  workable  system 
exists  whereby  differential  equations  may  be  simply  converted  into 
block  diagrams  and  vice  versa.    The  system  is  workable  for  various 
flight  conditions.    It  can  also  handle  disturbances  up  to  one  tenth 
radian  in  most  cases  but  is  restricted  to  linear  approximations. 

The  system  is  compatible  to  any  multi-loop  control  problem 
and  is  not  restricted  to  airframe  control.    A  standardized  form  now 
exists  upon  which  multi-loop  compensation  theory  research  may  be 
conducted  using  either  basic  hardware  or  differential  equations  as 
the  basic  system. 

The  graduate  work,  for  which  this  thesis  is  a  partial  require- 
ment, was  performed  while  the  author  was  attending  the  U.S.  Naval 
Postgraduate  School  during  assignment  to  the  Aeronautical  Engineering 
curriculum  (Avionics). 
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1 .  Introduction 

The  types  of  motion  that  result  from  a  disturbance  in  some 
equilibrium  flight  condition  and  the  transient  motion  of  the  aircraft 
in  response  to  control  movements  are  analyzed  in  the  study  of  the 
dynamic  characteristics  of  an  airplane.    Over  the  course  of  years  the 
effects  of  dynamic  stability  in  aircraft  have  been  known  to  designers  „ 
The  mathematical  computations  were  discussed  by  early  aeronautics 
pioneers  which  included  Lanchester,  Bryant  and  Glauert.    In  spite  of  the 
fact  that  dynamic  stability  effects  were  known,  most  people  paid  little 
concern  to  this  field.    As  design  groups  solved  other  problems  of  air- 
craft design,  the  dynamic  stability  effects  were  automatically  solved  as 
a  secondary  result.    However,  since  World  War  II,  aircraft  development 
has  pushed  into  the  field  of  very  high  performance  aircraft.    In  many 
cases  the  dynamic  stability  problem  doesn't  solve  itself  as  happened 
previously.    Furthermore,  the  fire  control  problem  has  entered  the  scene 
which  has  generated  a  control  problem  of  no  small  magnitude. 

This  work  will  consider  the  aircraft  in  six  degrees  of  freedom  and 
explore  paths  of  compensation  and  solution.    It.  is  the  prime  intention  to 
explore  the  work  of  Chu^'  in  multi  loop  servo  systems  and  determine  the 
value  of  it  as  applied  to  the  aircraft  stability. 

The  first  section  will  set  down  a  general  definition  of  terms  and 
the  basic  equations  of  motion  as  related  to  aircraft.    It  will  further  set 
down  the  basic  assumptions  upon  which  the  equations  are  formed  „ 
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The  equations  motion  for  the  airframe  which  are  completely 
derived  in  Appendix  I  are: 

<?  *  -  R7)m  =  X  -  m     sin  0  "-V 

s   (V   «.  pjj  .   pw)m  3  y  +  mg  00s  0  sin    O 

*   PV  -  QU)m  =  2  +  mg  cos  ©  cos    ^  (7 

L     -  AP  +  &K1  +   (G  -  P)H?   +  1  '  "  '   -  R2)   -  E(PQ  +  R) 

R  -  cK3:  -  bK0R  (T-U) 

M     ~  BC   +  bKo   +   (A  -  C)K    +  E(P9  -  R°)  -  FOR  +  DPQ 

♦  -K^R  -  cK3P  (1-5) 

=  CR  +  cKg  +   (B  -  A)P°  +  F(Q2  ~  P2)   *  E(QR  -  f) 

-  DI      *  '"K2P  «  ?.KX:  (1=6) 

He  +  Fe  -  h'i     s  Me*cz  -  ?ex  <RP  -  °->  <T-7> 

*  Fr  a  Jr  ?  s  mr1racy  =  <R  +  pM?rx  -  fR;  "p)rz  C1-" ' 

?  Hg   ■«■  ?a  -   Ia^"     +  2   Pay   (RC    +  P)  (1-9) 

The  ar  r   and  linear  velocity  of  a   free  body  in  space  ar^ 

p     ^  p~  if  sin  e  -10) 

-  e  cop  0    +    V  sin  0cos6 

P.     -  <f^os  ?5  co?  e-    ©sin  ^  (T~ 

dx*  -  ?  O  c^s  ^  +  V   (sin  0  sin  &  cos  ^  -  cop  £5  sirs  <f  ) 

dt 

+  W  (cos  0  sin  ©cor  ^    +  sin  ^  sin  tf5  )  (l-i^ 

dy*  -  TT  cos  0  sir  ^  +  V  (sin  C5  sin  O  sin  tf  +  cos  0cog^!) 
dt ' 

(r  >s  0  sin  O  sin  fy  -  sin  0  cos  ^  )  ( J=lli) 

It*  a  -TT  s?n  ^    +  V  sin  $  cos  B  +  W  cos0  cos  O  (i-l!?) 

dt" 


Consideration  of  the  aerodynamic  forces  must  be  resolved  into 
the  effects  they  produce  on  various  parts  of  the  body  in  some  flight 
condition  and  also  the  resultant  moment  effects  they  produce. 

The  first  or  necessary  condition  is  to  establish  a  reference 
condition  and  proceed  from  there.    A  reference  condition  is  some  flight 
condition  where  equilibrium  exists  in  all  equations  and  the  summation 
of  forces  is  zero.    This  is  a  steady  state  condition. 

The  set  of  axis  upon  which  this  rigid  body  is  oriented  is  commonly 
called  the  Eulerian  Axis  system.    General  aerodynamics  usage  modifies 
the  z  axis  so  the  positive  direction  points  downward. 
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The  first  rotation  moves  about  the  Z  axis  in  a  direction    '-f    and  is  known 


as  yaw.    The  next  is  the  movement  about  the  Y  axis  which  is  term     © 
and  known  as  pitch.    The  third  motion  about  an  axis  is  roll  about  the  X 
axis  and  is  maneuvered  in  terms  of   0    .    Caution  must  be  used  in  de- 
termining these  positions  because  the  order  in  using  this  system  is  not 
commutative . 

The  derivation  for  the  yaw  rate  R ,  the  roll  rate  P  and  the  pitch 
rate  Q  are  shown  in  detail  in  Appendix  Two.    They  are: 

P     =    JZS   -  ^sin  e  (1-16) 

Q     =    ecos  ^+   Cf'sin  ^  cos   O  (1-17) 

R     =    ^'cos^cose  -  &  sin    <ft  (1-18) 

and  where  small  angles  are  involved 

P     =    <fi  (1-19) 

Q    =      e  (1-20) 

R    ~     V  (1-21) 


-22) 
-23) 
-24) 


Similarly,  Etkin  shows: 

©    =  Q  cos  p  -  R  sin  ^5 

ft  =  P  +  Q  sin  £$tan  S  +  R  cos  <$  tan  3 

<jJ    =  (Q  sin^  +  R  cos   <fi  )  sec     Q 
The  general  equations  of  motion  have  to  be  applied  to  the  aircraft 
motion.    The  effects  of  axis  movement  upon  the  weight  force  of  the 
aircraft  will  be  shown  as  well  as  the  reference  state  condition.    It 
must  be  remembered  that  the  gravity  force  is  always  oriented  toward 
the  center  of  the  earth.    The  reference  body  axis  may  or  may  not  be 


oriented  in  this  direction 


From  the  discussion  of  the  Eulerian  angle  concept,  it  is 
determined  that  a  reasonable  reference  condition  would  be 


F    =  O  =  X     -  mg  sin  © 
x  o  ° 


Fy  =  Y    -  W  cos  e       sin  0  (I- 

o  o  r    o 

Fz  =  0  =    Zc  cos    &a    cos    0O  (I- 

This  allows  for  a  steady  state  condition  where  the  aircraft  is  not 


accelerating  but  neither  is  it  necessarily  in  straight  and  level  flight,, 
From  the  previous  illustration  of  Eulerian  angles  coupled  with  the  deri- 
vations in  Appendix  II,,  the  gravity  forces  can  be  determined  as  shown 


below. 
ccs  G cos  Q-* 


cosOsin  <f 


-sin<S 


cos  <P  Pin  ©sin  c6  cos  ^  cos  <fi  cos  O sin* 

-sin^cos    <p  +  sin1^  sinO'sin  0 

costf  sin  ©cos  <jf>  sin^sin©cos  (p        cos©cos0 

-  sin 4*  sin  <fi  -cos  ^  sin  d 


my, 


y 


myr 


J 


Discussion  and  derivation  of  this  general  cosine  matrix  are 
found  in  several  sources  (2,3) .    The  overall  result  of  the  equations  oi 
motion  is  shown  below. 


- 


O  Q        r^+  mp    (cQS0Q   sin/$0) 

S         V"    -  '     ^''         s^  )(cos©?in^) 

ir>©    Kcostp  sin©  sin  ^  -  sin  <f  cos(Z$) 

+  trig   (cos©     sin 0O) (cosy cos  C25  +  si^ty  sin  £>'=;>->^,/ 
+  mg   (cosOc  cos^    )(cos©sin^) 

m   (W  +  FT"  -  ' TT)   •«  '"    -  m?r   (sin^Q)(cos  </*sin  Scos0  -  siny'sind&) 

+  ■  ^©0  sin^0)(sin^sin©cos^5  -  cos^sin<Zi) 

+  m£   (cos^0  cos ^Q) (cof, Ocos ai)  -30) 

L  -  A?  +  aKx  +     H   (C-B)   -  r-   (PQ  +  P^   +  c  K3Q  -  b*  R  -31) 

M  =  BQ  +  bK2  +  PQ  (/-C)   +  E  (P2  -  -?)   +  -K]R  -  cKP 
3E   *  e  K^    *     '     (BW  )   +  E   (QR  -  P)  +  bK«P  -  -  Ki '"• 


The  above  equations  are  complete  except  for  the  external  forces 
X,  Y,  and  Z  and  the  external  moments  L,  M,  and  N.    These  forces  and 
moments  are  composed  of  aerodynamic  and  propulsion  effects  and  the 
moments  due  to  control  surfaces.    These  equations  are  very  general  in 
nature  and  contain  the  following  assumptions  for  simplification  frorr   the 
most  general  and  complex  case.    These  assumptions  are: 

1.  The  earth  is  assumed  to  be  fixed  in  space  and  the  Earth's 

atmosphere  is  assumed  to  be  fixed  with  respect  to  the  earth, 

2.  The  airframe  is  assumed  to  be  a  rigid  body. 


3.  The  mass  of  the  body  is  assumed  constant  during  the  time 

of  the  problem. 

4.  For  an  aircraft  or  missile,  the  XZ  plane  is  assumed  to  have 

mirror  symmetry . 
The  above  equations  then  cover  the  six  degrees  of  freedom  found  in  an 
aircraft  or  missile.     Coupled  with  the  other  auxiliary  equations  covering 
roll,  pitch  and  yaw  effects  along  with  control  motions,  solution  is  pro- 
bably possible.    To  date,  this  has  been  accomplished  either  by  a  machine 
using  specific  coefficients  for  various  terms  or  several  simplifications  „ 
Discussion  of  this  will  be  conducted  in  succeeding  sections. 


2.  Consideration  of  External  Forces  and  Moments  on  a  Body  in 

Motion 

In  the  previous  section  the  equations  of  motion  were  discussed. 
The  aerodynamic  thrust  and  control  effects  were  simply  grouped  in  Xt 
Y,  Z,  L,  M  and  N,    It  is  the  aim  of  this  chapter  to  discuss  these  in 
detail  with  respect  to  the  effect  on  the  dynamic  control  stability  problem 

Consider  the  term  X  when  the  aircraft  is  in  a  disturbed  condition. 

The  reference  condition  is  X    and  thus  X  =  X    +  AX  for  some  disturbed 

o  o 

condition.    The  question  is  then  to  determine  the    A  X  components . 

The  force  and  moment  variables  can  be  expressed  in  a  Taylor's 
Series    . 


F  =  Fo  +  /^?W+/§M  *1  —  +  (9  f)/q 


o2?  A2 


3fi 


9» 


For  practical  reasons  the  effects  of  second  and  higher  order  terms  are 
omitted  from  further  considerations  as  will  be  discussed  later.    Con- 
sidering the  motion  of  the  aircraft  itself  the  X  forces  would  be 


X  «  Xo 


m*ffi*mww^  (ii"2> 


Similarly , 


•     • 


wrwmmw  m11'3 


For  convenience  all  terms  will  be  expressed  in  letter  subscript  form. 
It  can  readily  be  seen  then  that  considering  X,  Y,  Z,  L,  M  and  N  and 
the  effects  of  U ,  V,  W,  P,  Q,  and  R  plus  the  acceleration  terms  U, 
V,  W,  P,  Q,  and  R  there  are  72  derivative  terms  required.    This  does 
not  include  the  control  effects  which  can  add  another  36  derivatives. 
Needless  to  say,  considerable  effort  has  been  expended  in  justifying 
simplifying  assumptions.    In  addition  to  this  ,~  obtaining  reliable  co- 
efficients for  various  derivative  terms     is  not  always  possible. 

It  is  not  the  intention  of  this  work  to  cover  the  computational 
methods  for  these  various  derivatives.    The  scope  is  too  large.    A 
search  of  all  authoritative  books  [3,  4,5,  6  Jon  the  subject  yields  a 
table  of  terms  which  can  be  reasonably  determined.    This  appears 
Table  I.    These  considerations  are  mentioned  to  bring  out  the  awareness 
of  the  different  causes  and  effects  of  external  forces.    There  are  con- 
siderable gaps  in  the  table.    Determination  of  these  to  complete  the 
table  would  be  a  topic  of  a  complete  study  in  itself.    Much  work  remains 
to  be  done  in  this  field. 

Furthermore,  some  of  the  methods  of  computing  the  known 
stability  derivatives  are  questionable.    Some  computations  from  theory 
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prove  quite  accurate.    Others  are  very  poor  and  empirical  methods  or 
model  testing  are  used.    Even  these  fail  at  times  and  "educated  gues   e: 
from  previous  airplanes  are  used  as  a  basis  of  determining  stability 
derivatives  for  a  new  design.    Tail  effects  are  the  most  notable  in  thi 
last  category. 

In  all  aerodynamics  problems,  considerable  effort  is  made  to 
non-dimensionalize  the  equations.    The  advantages  of  doing  this  are 
of  considerable  merit.    For  the  most  part,  it  is  rather  simply  done. 
It  has  been  found  in  both  theory  and  actual  practice  that 

F=4/3V2SC  II-4 

Z.  F 

Where  S  is  a  reference  area,  usually  the  wing  area,  V  is  the  absolute 
reference  velocity  and  Cp  is  some  non  dimensional  coefficient.    C 
varies  with  Mach  Number  and  Reynolds  Number  depending  on  the  force 
considered.-  (lift,  drag,  etc.)  r     is  the  density  of  the  I  luid  through 
which  the  body  is  passing.    Thus  a  term  such  as  Zw  has  the  following 
relationship: 

zw1  -  Jl         ^2s  /  a  °'7  V  */°v23cZw  11=5 

Here  again  V  and  S  are  the  velocity  and  wing  area  terms  respectively. 
The  same  concept  applies  to  the  moment  terms  except  that  a  moment 
arm  term  must  also  appear  in  the  equation.    A  generalized  example  is: 

M  -  />  *  !I"6 

In  this  example  M  is  some  moment  and  C      is  a  needed  dirnensionle 

M 

coefficient  to  satisfy  the  equation.     P  ,  V  and  S  are  the  same  as 
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mentioned  above  while^is  the  moment  arm. 


Likewise  in  a  stability  moment  the  derivative  can  be 


form . 


1 


(f? 


»  f~T'31^ 


Another  external  force  to  be  considered  is  that  of  thrust. 
This  becomes  an    important  factor  in  most  flight  conditions.    The 
location  of  the  thrust  axis  is  seldom  coincident  with  any  reference 
axis  in  the  body. 
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Assuming  the  thrust  line  lies  in  the  xz  plane,  then  the  thrusl 

forces  can  be  resolved  in  this  generalized  manner,, 

X     =  T  cos  6  II-8 

T  T 


Z     =  T  sine^  II-9 

T  T 


MT  =  T7  11-10 

For  a  steady  state  condition  the  equations  become 

X0    =  To  cos  O  T  II-ll 

T  ■*• 


ZQ  ,  =  To  sinB  T  11  =  12 


IVI        =    T0lT  11-13 

Om  °Z 


while  in  a  disturbed  condition  where  the  body  axis  remains  fixed  in 

the  mass , 

Xm  =  T.  cos  e  „,  11-14 

T  1  T 


ZT  -  Tj  sin  e  T  11=15 


MT  =    Tlz  11= 1G 

and  T.  =  To  +    AT  11-17 

1 

Assuming  the  air  density  remains  constant  so  the  engine  thrust  will 
not  be  affected 


\d&      J 


11-18 
or 
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These  are  the  considerations  for  the  external  forces  on  a  b< 
in  motion.    While  the  dimensionless  coefficient  has  been  used  exten- 
sively by  the  aeronautical  engineering  field,  it  will  not  be  used  furl 
in  this  work.    The  primary  mission  of  this  work  is  to  investigate  the 
control  response  and  paths  for  possible  compensation. 

The  generalized  equations  of  motion  then  appear: 

a 

I  -  W)  +  mg  sin60   (cos©cos4      -  osOQs±r\^ 

S0    )(s*r\&)   -  X     +  X1   u  +  X1   w  +■  X**w  +   •'  '  ,~n 

r  o   v  o  u  w  <■:  o  \ 

\K      '  •  coseT  +  T"  C°seTu  +  Tcfrpm  cos^l£r3 

v    i  rag  sineo(cos^sin£sin^~  sin^cos^)  »    i 
(cos<^cos^  tyi      e>       0s   -  6  cos^  )(c-,^B     -  -7: 

-v  +  Yfr         -      r  +  Y ■*     c    +  Y»     V  "-^ 

>'ipeo(cc,r4>s'->ecos^  +  sin^sin^T) 
a     :  -  ^5;iV-:nff-3in©cos^  -  ros^slnS)  -  mg(i    >    6  ^  ' 

©     s(2$  '    i        •  w  +     »•  '  a  +  *.»  t?    +  z6o"yf 

K-,  '  +  ck       ■-  '•  K~R  =  L     +  L8   v       L«*v 

c  O  V  * 


^  ftF+L'^  '  s 


K  ;<2  -  R2)   +  aK  R  -cKP  +  M*  u 

1  *  u 


'  *"  •  ' )         (QR-P)  h      K„ p  -  s  K-. 

1 


*§    ^    £ 
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Thus,  it  has  been  shown  that  the  equations  of  motion  have    the 
above  external  forces  which  will  affect  the  rigid  mass  in  motion 0    The 
assumptions  made  are: 

1.  The  earth  is  assumed  to  be  fixed  in  spaee  and  the  Earth0 
atmosphere  is  assumed  to  be  fixed  with  respect  to  the 
Earth . 

2.  The  mass  is  a  rigid  body. 

3.  The  mass  of  the  body  is  assumed  constant  during  the  time 
of  the  problem. 

4.  For  an  aircraft  or  missile,  the  xz  plane  is  assumed  to  have 
mirror  symmetry. 

These  assumptions  were  made  in  the  first  section,    in  addition  the 
following  assumptions  and  simplifications  are  made  in  considering  the 
aerodynamic  forces . 

1 .  The  flow  is  assumed  to  be  quasi-steady. 

2 .  The  effects  of  auxiliary  aerodynamic  equipment  such  as 
flaps,  speed  brakes,  slats  and  spoilers  are  neglected „ 

3„      The  thrust  can  be  affected  only  by  change  of  speed  and/or 
revolutions  per  minute.    No  change  is  considered  for  change 
in  fuel  mixture,  manifold  pressure,  side  slip  or  propeller 
fin  effect. 
These  assumptions  were  made  for  several  reasons.    In  the  first  and 
third  cases ,  absence  of  adequate  information  on  many  of  the  coefficients 
makes  ignoring  or  setting  them  equal  to  zero  mandatory.    In  ail  probal 
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these  coefficients  are  very  small.    Assumption  two  covers  speci  A 
cases.    While  these  items  may  have  considerable  effect  when  used;, 
they  are  eliminated  from  the  general  case  at  this  time  because  an 
craft  probably  uses  them  less  than  two  percent  of  its  flying  time*    For 
the  special  cases,  these  effects  may  be  quickly  added. 

The  work  in  these  first  two  sections  is  developed  primarily  for 
background  purposes.    There  is  nothing  new  or  startling  in  this.    For 
further  examination  of  this  subject,  the  reader  is  referred  to  the  woi  • 
in  Ref .  7  and  8  which  closely  parallel  the  development  in  this  work. 

The  aerodynamic  force  and  moment  terms  carry  a  primed  notation 
for  ease  in  publication.    In  later  sections  these  terms  will  be  divided 
by  mass  or  moment  of  inertia.    These  terms  after  division  will  not 
carry  the  prime  C  )  notation .    There  are  far  more  of  the  latter  terms 
used.    Thus,  it  is  for  ease  of  publication  only  and  does  not  affect 
theory  in  any  way. 
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3.  The  Airframe  Equations  in  a  Standard  Form 

This  section  will  take  the  general  equations  developed  in  I 
previous  chapters  and  simplify  them.    The  simplification  will  exten 
to  an  aircraft  in  straight  and  level  flight  and  subject  to  small  per- 
bations .    This  is  the  standard  method  of  dealing  with  the  problem  as 
found  in  present  day  texts.    However,  these  equations  will  be  use 
in  block  diagram  form  and  in  determinental  arrays  .as  developed  in 
Chu°s  work.  ^--)    The  results  will  be  compared  with  present  day 
classical  methods  to  determine  validity  of  this  determinantal  concept „ 

The  first  simplification  is  the  concept  of  the  lateral  system 
being  independent  of  the  longitudinal  system.    Therefore,  the  six 
equations  of  motion  are  divided  into  the  two  groups.    The  longi- 
tudinal group  includes  the  directional  motion  in  the  X  and  Z    directions 
and  rotary  motion  about  Y  axis.    The  lateral  group  consists  of  direc- 
tional motion  along  the  Y  axis  and  rotary  motion  about  the  X  and    Z 
axis .    The  axis  system  referred  to  here  is  the  body  axis  system  of 
the  airplane. 

Secondly,  the  aircraft  is  initially  in  a  steady  state  condition 
This  is  the  reference  state  as  mentioned  previously.  When  the  dis- 
turbance terms  are  all  zero  Equation  11-19  reduces  to; 

Up  +  QoWo  -  RoVo)   +  v.  B  0  ■  Xo  +  T     cos  e  T  III-l 

However,  it  was  stated  that  this  situation  was  restricted  to  longi- 
tudinal or  lateral  systems.  Thus  for  motion  along  the  X  body  axis , 
a  yaw  motion  Ro  and  a  side  motion  Vo  would  be  zero.    Further  it  v 
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stated  that  the  aircraft  was  in  a  steady  state  reference  condition  „ 

a 

Therefore  y  would  be  zero.    Thus  equation  III  —  1  reduces  to 

m  (QoWo)   *  rog  sin  Sn  *  X'c   =  To  res  6    T  jjj. 

Similarly,  it  may  be  observed  that  equation  11=20  for  the  steady 

state  reduces  to 

Yo  =  0  III- 3 

Likewise  equation  II™ 21  reduces  to 

m  (QoVo)   --»  Z»     +  mg  cos  ©      +  To  s:n  6  T  -  0  HI-4 

o 

The  moment  equations  all  reduce  to  zero  provided  the  rotor  effects 
are  neglected  or  considered  zero. 

Lo  =  Mo  -  No  -  0  III-5 

Recalling  that  only  small  disturbances  were  to  be  considered 
and  these  disturbances    were  to  be  in  only  the  symmetrical  or  longi- 
tudinal system  for  the  longitudinal  -group  the  derivative  due  to  asy- 
mmetric variables  of  V,  P,  R,  £  and    c;    are  zero.    In  a  similar  manner, 
when  the  asymmetric  system  effects  are  considered  the  symmetric  d 
turbances  are  ignored.    This  leaves  the  variables  V,  W,  Q ,  and  Y[ 
out  of  the  a  symmetric  equations.    Finally,  because  of  the  small  dis- 
turbances, the  change  angle  Gwill  be  small  so  the  approximation 
cos  0  =*  1  and  sin  0   =    J0       is  valid. 

The  result  of  the  previous  discussion  is  subtracting  equations 

III— 2  from  11-19,  eliminating  the  asymmetric  changes  and  supply 

the  small  disturbance  considerations.    The  result  is: 

m  (u  +  ">*o  E*>%^4  X'  +  n»fi  e  ~o?  e  0  *  :ru»u  +  x^w  +  x^»w 

III 
x'n  *  ::^,lTt  +  T-1  coseTu  +  T^»RIMcoseT»<$i 
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For  the  z  direction  taking  IXI-4  from  IX -21  and  using  the  same 

assumptions: 

V     -  \  ,;.■■  ITi 

For  the  moment  in  the  y  axis  „  the  following  applies: 

11  "  I  III-8 

*  T  y    I    &     - 
To  complete  the  case  for  the  symmetrical  group,  the  action  of  the 

elevators  must  be  included.    From  appendix  I  the  elevator  hinge 

moment  effects  are  reduced  to: 

"  »         i  '*  t  i 

Tp-n  =  tr  !  r?+  H  '.Tt  v&-  *  FeiiU  *  HecO  +   il^r  III- 9 

In  most  cases  where  moment  of  inertia  terms  are  handled  It 
is  usually  expedient  to  use  an  axis  system  which  is  the  principal 
inertia  axis.    This  means  that  the  product  of  inertia  terms  are  zero. 
For  aerodynamic  motion  problems ,  this  applies  for  some  cases  „  ^  " 
This  case  is  not  one  of  them.    What  is  gained  in  simplification  of  the 
rotary  motion  equations  by  use  of  principal  axis  is  lost  in  attemptin  \ 
to  handle  the  force  equations  with  the  wind  velocity  consideration 
in  sine  and  cosine  terms .    By  orienting  the  X  axis  into  the  reference 

wind  velocity  the  result  is: 

7  2  =  ' 

III- 10 

where       U  =  U    +  u . 

Therefore,  in  the  steady  state  reference  condition V     =  U     and 

c         o 

V    =  W    =  0.    With  the  assumption  that  u,  v  and  w  were  very  small, 
their  products  and  squares  can  be  neglected. 
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and  since  U    L  is  much  greater  than  2  Unu,  then  V  "=U    ^u  and  the 

terns  are  interchangeable  in  this  section  . 

-\    previously  stated,  the  reference  situation  is  a  steady  si    t€ 

motion  and  includes  first  order  effects  only.    The  orientation  of  the 

axis  U     is  a  finite  value  while  \/    =  W    =  0  „    Furthermore ,  the  author 
o  o         o 

is  hard  pressed  to  visualize  an  aircraft  or  missile  tumbling  along  v. 

a  Q~  motion  although  P     exists  in  spinning  missiles  and  bullets  „ 
u  o 

Further  X     s  X    =  Z  *  =  X  >?  -X-vi    =  T_  nniv ,  =  0  as  an  assumption  0 
w        q         w  i  S 

These  terms  are  usually  very  small  „ 

The  "standard"  determinant  form  as  stated  in  Chapter  3  of 
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Chu1  is: 

1   -  FF! 

*n 

PF2 

PF21 

-  PF, 

1   +  PF? 

FF22 

■    ] 

-  PF2 

FF3       PF30 
1  *  FF3 


where  the  PF  indicates  a  performance  function .    The  system,  may  h 
/ye,  nodes.    This  example  contains  only  three  nodes.    The  "standard" 
block  diagram  for  this  form  is: 
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are  3-3        TANDARD  *LOCK  DIGRAM  DRAW  FROM  3x3  DETERMINANT 
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The  nodes  are  labeled  in  Roman  Numerals,,    The  basic  rules 
of  nodal  signals  allow  several  inputs,  but  only  one  output.    T      - 
output  fro:::  the  node  goes  directly  to  a  performance  function  with 
single  subscript  number.    Signals  may  be  picked  off  fro":  the  node 
output  and  fed  forward  to  another  node.    An  example  of  this  is  contained 
in  the  signal  going  through  PF.    .    In  this  case  the  signal  is  also  mull 
plied  by  the  function  PF-»o  and  then  fed  to  node  III.    Feedback  tak« 
place  in  a  similar  manner.    The  signal  output  from  the  first  system  PF 

JL 

is  picked  off  and  fed  through  PF-,,  back  to  node  I.    Similarly  the  other 
nodes  function  in  the  same  manner. 

For  feedback  blocks ,  the  first  subscript  number  denotes  the  out- 
put signal  source  for  the  feed  back  and  the  second  subscript  number 
the  input  node  to  which  the  feedback  signal  is  going.    The  notable 
points  about  the  determinant  are: 

1 .  The  main  diagonal  always  has  terms  of  one  plus  the 
direct  path  performance  function  times  the  self-loop 
feedback  performance  function . 

(1  +  PF1PFn) 

2 .  All  terms  above  and  to  the  right  of  the  main  diagonal 
are  feedback  terms . 

3.  All  terms  below  and  to  the  left  of  the  main  diagonal  of 
the  "standard"  determinant  form  are  of  a  feed  forward 
nature.    All  of  these  feed  forward  terms  are  negative. 
The  term  immediately  below  the  main  diagonal  term  i 
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direct  path  performance  function.    Thus,  where  the  second 

column  main  diagonal  spot  has  I  +  PF  PF„0if  the  direct  [ 

2      *^ 

term  PF?  will  appear  in  the  location  immediately  below 

-PF2. 
These  rules  are  strictly  mechanical.    For  a  rigorous  derivation  of  th 
system  the  reader  is  referred  directly  to  Chu°s  work. 

In  Chapter  9  of  Etkin  *  ' ,  there  exists  a  two  degree  of  free 
situation  where  ©and  W  vary  and  the  speed  in  the  forward  direction  is 
constant.    This  is  the  short  period  approximation.    The  equations  us 
the  LaPlace  transforms  and  nondimen  ional  coefficients  are; 

)  aC   -    ft>  .      +  (  &  •       Y?      a   0  HI 


rj  )e        -t.17  ni-12 


5-      .©-  ?hn,7l*A  m": 


This  author  also  grouped  his  acceleration  terms  such  that  G^ 


included  (  and  G  i  .    Similar  situations  exist  for  Gv^     ,  Gn  #> 

and  G  .    Note  that         <f  and  p       P    for  this  consideration .    Froi 

these  equations  a  block  diagram  is  evolved. 
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Further,  Etkin  derives: 

n  >Zok)(iBS2-G       )  -Gm  I 


®n 


111=14 


G7   )s 


J*>{ 


A0i,°z,'V?-S)^^,^)8 


III 


n     (2 


^: 


z^ifi-32-^-0^^*0 


Not  it  is  intended  to  use  Chu's  work  involving  the  use  oi 
"standard1"  block  diagrams  and  "standard"  determinant  form  to  prove 
that  the  "standard"  system  applies  for  the  airframe  as  well 
multiloop  ser/o  system. 
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The  Z,  M  and  "»j   equations  of  motion  with  the  applied  assump- 


tions and  $     =0  are: 
o 


v   «    'T. 


"l7  ^ 


III- If 


-y[  -  -n^l 


111=17 


rf-  io^  ~  \z\\ 


10    J  -  10  *;  -  10  „q 
w 


-    A  Fe 


III 


Recalling  that  u>  =  Uo   tan  A  o^and  for  small  angles 
•a)  ~  U0  A  <*,   and       ©  '=*  q  8 

The  equations  in  LaPlace  notation  then  become: 
(s2  =■  (sl0-»J  *  10iO)7£-  (sl0w-  *  10w)  U0A<*«  (slOq)e-AFe  III-l: 

(s2  -  sl\,)G  -  (sM>J  +  Mt?)7/  -  (sM^*  1^)  tjoAK        -  0      III   : 


(s  -  zw)  u0A<x  -  (slJ0  *  c"o)e  -  (s2  *i  *  7  l]\m  ° 


III 


These  equations  in  standard  determinant  form  now  become: 


\ 

A  < 

i 

3   -  (sIOt?   ♦  10n) 

3 

• 

-(slOoC    *  10H  ) 

S1" 

-  s30o 

s 

-(.zi  »z,) 

1  -  L* 

-  s(Uo  + 

5 

s 

s 

• 

-  (sM*<  +  M«0 

2 

s 

1  -  sMg 
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Ordinarily,  it  is  not  necessary  to  divide  by  s  or  s    .    How- 
ever, it  is  done  here  to  demonstrate  more  closely  the  standar 

To  take  the  differential  equations  and  arrange  them,  by  rowi  and 
columns  to  see  the  actual  physical  system  is  to  no  avail  „    The  an 
ment  doe3  not  exist  in  the  above  form„    It  will  be  shown  how  to  achieve 
this  form  later. 

Consider  the  "standard  form"  array  as  set  forth  above.,    The 
block  diagram  for  this  in  Chu's  "standard  form"  is: 
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;ure  3-3     RE-STH.T  OF  BLOCK  DJ'  NATIONS 

?s20  and  21 
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The  above  block  diagram  then  is  a  three  node  system  but  the  o  itp  i1  . 
from  the  first  and  second  main  path  blocks  are  nothing  that  is  n 
recognized  in  an  airframe.    Because  of  the  fact  that  the  main  block 
following  the  third  node  is  1 ,  the  output  here  is   @    „ 

A  few  simple  block  diagram  manipulations  clears  the  situation 
and  shows  true  signal  outputs  that  are  physically  realized  in  the  air- 
frame.   Block  diagram  manipulations  allow  the  changing  of  a  pick  off 
point  provided  the  transfer  function  is  corrected „    Further,  two  bloc 
can  be  combined  into  one  equivalent  transfer  function. 


-(si  AH    v, 

n  n 


which 

changes 
to 


[AGRA! 
Adding  a  dummy  node  this  now  changes  to: 
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The  result  of  using  these  two  maneuvers  rearranges  the  block  d 
to  the  following" 


OF  F  SURE  3-3,   US] 

OF 


- 


Maneuvering  the  above  diagram  location ,  the  block  diagram  1    *  e     ; !  -. 
form  of  Fig.  3-7,  which  coincides  with  Fig.  3-2. 
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FRKEna 

re  1-6  to  compare  with 
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It  is  noted  that  this  block  diagram  approaches  the  form  set 
down  in  Etkin  except  for  nondimensionalizing  the  terms.    The  dashed 
area  in  the  diagram  is  what  Etkin  labels  Gq^        and  G  o(ti 
.-The  determinant  from  the  above  dashed  portion  of  the  block  di  igra 
is: 
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Evaluating  for  the  transfer  functions  -&   and  ^^- ,  the  results 


will  be  in  the  form     ^rjk^iy 

(1) 
This  form  is  the  type  expressed  by  Chu       and  Thaler 

term,  which  is  the  whole  determinant,  reduces  to: 


The  A 
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^r  4  =  -   -  V  )(s5    -  sMq)   -  s(Uo  +  Zq)(sM^    +  Mo() 


III-2; 
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To  consider  the  &  oc    term  the  cofactor     A  y  of  this  determinant  is 


o3  t    ined  by  crossing  through  the  first  row  and  the  second  column 
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The  m;dn  feed  forward  transfer  function  between  the  input  signal  fro: 
node  2  to  the  output  is  +1  in  this  case.    This  is  also  true  of  node  4 
input  to  output  ,    Hence  G»2&<*     =  +1. 
The  transfer  function 


6  « 
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'e-7 


A 

^    *  Z^)(s2  -  sM  )   4  s(Uo  -  7q)(sMi  +  M^)     IH 


(s  -  ?oO(s?  -  sMq)   -  s(TJ0  *  ZqKsM£  +  Mc<) 


It  can  be  noted  that  this  equation  111-24  is  the  same  equation  as 
III—  1 5  given  by  Etkin  save  the  non-dimensionalizing  terms. 
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To  evaluate  the  transfer  function  &  i  7)   ,  the  same  method  is  us 
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&  '  *  ri  Vs    -  £<*)    +   (sMci    ♦  McC)(s*?>    <    Z-w   [II 


.-■^(•c   ,    ^      +  lSKo^  +  iic^  ; ,.?  •  -^   < 


i 


Therefore: 


(     ^ 


vS   .  Z*)(s2  -  -   s(lT     +  Zn)(sMoC  +  Mo<) 

I  °  - 


Here  again,  excepting  the  non-dimensionalizing,  the  result  is  the 
same.    Equation  III  — 26  is  equal  to  equation  III— 14  as  given  by  Etkin  . 

The  system  of  "standard"  forms  then  works  for  the  airfran  e 
differential  equations  after  some  manipulation.    It,  can  be  noted  that 
there  is  a  method  of  making  up  this  standard  form  without  the  necessity 
of  the  manipulation  which  will  yield  physically  realizable  output  sign 
When  the  block  diagram  of  Fig.  3-6  is  converted  into  a  determinant 
for  solution  of  the  transfer  function,  there  appears  a  series  of  alterna- 
ting rows  where  the  equation  -1+1  =  0  appears .    It  also  appears  t 
these  simple  equations  occur  in  a  sequence  where  the  plus  one  ten 
placed  in  the  2n!  row  and  column,    n  is  the  number  of  differential 
equations  in  the  problem.    The  minus  one  term  is  placed  in  the  2nl  row 
and  (2n!"l)  columns.    All  feedforward  terms  (those  below  and  to  the 
left  of  the  main  diagonal)  are  moved  one  column  to  the  right  and  rei 
in  the  r-  •*  e  row.    This  doubles  the  size  of  the  determinant.    Hov. 
there  are  several  zero  terms  introduced  and  the  actual  work  involved  in 
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evaluation  of  the  determinant  is  not  markedly  increased . 

The  next  item  is  to  look  at  the  lateral  motion  proble: 
slightly  more  complex  because  there  are  five  differential  equation 
c:      otion  instead  of  four  „ 
Equation  II— 2  with  simplifications  becomes; 

r)   -  flmz  ~c&Oo  -  :  Yp  -  YRR  .  Y.    J  8  0  111-27 

0*  "^  o "  \fi~  Y  pp  -  Y  -r  -  Y£  ^  s  °  in~28 

—^  _    _    __ 

The  moment  equation  11-22  becomes 
P-SH.  T,y5    -  Ly  -  Lpr  -  L.|     -  L^J  -  0  III=29 


and  equation  11-24  becomes 

R  +  EP  «  W   fi     -  NpP  -  NpR  -  N^  £     -  N*    £  -  n/  £    *  0 


III- 
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The  aileron  hinge  motion  equations  from  Appendix  I  are 
£  -  UpP  -  ll^r  -  lie  |    -  11  v   J    -  A  F* 

And  also  from  Appendix  I ,  the  rudder  hinge  motion  is 
\  -  1  V>    ~  12PP  -  12rr  -  1?£  £  -  1?^*  £    «  A  ?r  Hi 
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From  these  terms  a  determinantal  array  can  be  made.    It  she 
be  noted  that  this  time  a  set  of  equations  of  -1+1  =  0  will  be  insei 
in  the  system.    This  permits  direct  drawing  of  the  block  diagram  and 
omitting  the  need  of  manipulation.    Table  3-1  then  presents  this 
From  this  Fig,  3-8  is  drawn  directly. 
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Figure  3-8  BLOCK  DIAGRAM  OF  AIRFRAME 
LATERAL  STABILITY  SYSTEM  AS  DRAWN 
DIRECTLY  FROM  TABLE  3-1 

(Roman  Numerals  indicate  node 
riumber) 
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.  ie  resulting  block  diagram  very  clearly  parallels  that  o!  tl  e 
lateral  system  shown  by  Etkin  in  Chapter  9^.    The    '  stand  to 

achieve  this  block  diagram  differs. 

There  Is,  then,  the  question  of  why  does  this  differ  from  r  <.. 
original   'standard'"  form  set  forth  by  Chu. 

Chapter  Two  of  Chu**'  discusses  at  some  length  the  idea  of 
arrangement  of  varying  performance  functions  into  a  ce standard "  blo<  • 
diagram..    The  individual  pieces  of  hardware  are  at  hand  with  known 
individual  performance  functions.    The  orderly  arrangement  of  these 
permit  the  construction  of  a  "standard"  determinant  for  the  whole  syatei 
From  this ,  a  multiloop  servo  system  can  be  handled  in  generalised 
results.    In  the  case  of  the  airframe,  the  determinant  formed  by  the 
differential  equations  was  known  and  the  block  diagram  had  to  be 
formed c,    Admittedly  a  block  diagram  (Fig.  3-3)  was  made  up  but  the 
output  signals  are  not  the  physically  realizable  type  that  are  generally 
known  to  pilots  and  engineers.    The  same  performance  function  resu  ,1 
from  the  system.    Both  are  mathematically  correct. 

The  difference  lies  in  the  fact  that  Chu  is  taking  Individ.. 
piece,   of  hardware  where  the  functions  are  known.    The  blocks  are 
made  up  and  then  the  "standard"  determinant  form  is  made.    In  the  <    i   e 
of  the  airplane,  the  equations  of  motion  are  known  and  a  determinant  is 
immediately  formed.    Each  column  contains  the  performance  function 
multiplied  by  the  direct  function  signal  output  term.    Chu's  standard 
forr  did  ite  s  that  all  terms  above  the  main  diagonal  are  multiplied 
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the  node  output  or  input  to  the  direct  function.    Steps  must  be  t  tken 
then  to  change  this  in  the  differential  equations.    The  addition  of  the 
alternating  -1  +  I  =  0  equations  solves  this  very  nicely    ind  in  no  v.    . 
alters  the  value  of  the  equation, 

A  series  of  evaluations  of  determinant  forms  is  given  m  Appen 
III.    The  results  show  that  the  addition  of  these  terms  in  no  w  ly  cl    n  ... 
the  final  value  of  a  determinant  or  a  cofactor  value™ 

Referring  to  Fig.  3-6,  there  are  other  block  forms  that  may  ch   i 
the  values  within  the  blocks,  but  will  not  alter  the  final  outcome.    Con- 
sider the  first  node  in  Fig.  3-6  as  an  example.    The  four  forr  is  expres   e 
in  Fig.  3-3  are  equivalent.    The  original  form  is  shown  in  Fig.  3 
The  self-loop  feedback  path  is  combined  with  the  direct  path  in  Fig, 
3-9 (b).    The  s     term  is  eliminated  and  the  performance  functions  are  then 
in  the  for-    shown  in  Fig.  3-9 (c) .    A  further  change  can  be  made  by 
dividing  the  nodal  input  terms  by  the  denominator  of  the  direr  t  ;  v  h 
feed-forward  function.    The  result  is  then  in  the  form  shown  in  Fig„ 
3-9{d)c    The  performance  functions  in  each  block  differ  somewhat,,  but 
the  final  result  remains  the  same.    These  manipulations  allow  the  use] 
more  freedom  in  the  forms  of  the  performance  functions.     It  illcv 
change-  in  form  to  fit  the  needs  of  the  particular  probler  . 

I  he  pick-off  point  for  the  feed  forward  terms  does  differ  in 
requirements.    It  is  seen  in  Figs.  3-9 (a)  and  (d)  th  i1  the  pick-oft  c  an 
be  either  off  the  output  signal  of  node  I  or  II .    The  final  result  will  not 
change,    However,  this  is  not  the  case  in  Figs.   3-9  (b)  and  (c) .    The 
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oint  must  be  off  the  output  signal,  of  node  II =    The  fin  .. 
re      It  will  differ  if  the  output  of  node  I  was  tapped  in  the  latter  c 


3-9(a) 


Figure  3-  : 
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For  convertin  ;    »  set  of  linear  differential  equations  into 
standard  block  diagram  form,  the  following  rules  can  be  applie 

1.  Between  each  differential  equation,  an  equation  -1  + 

ill  be  inserted.    An  additional  equation  -1  +  1  =  0 
be  added  at  the  bottom  of  the  set  so  the  number  o: 
tional  equations  shall  equal  the  number  of  different 
equations  originally  in  the  problem . 

2.  The  arrangement  of  the  values  of  these  added  equations  in 
the  determinant  shall  be  such  that  the  positive  one  value 
is  always  in  the  main  diagonal  block.    The  minus  one  term 
shall  always  be  placed  in  the  block  immediately  to  the  left 
of  the  main  diagonal  term  . 

3o     All  other  blocks  in  the  rows  of  these  added  equations  shall 

be  of  zero  value . 
4  „     All  finite  performance  functions  of  the  original  different- 
equations  which  are  located  below  the  main  diagonal  sh  il3 
be  moved  one  place  to  the  right.    The  vacated  block  5  si     1 
be  replaced  with  a  zero  value.    This  last  rule  is  option  , 
to  the  user.    It  does  save  the  possibility  of  error  which 
could  appear  because  of  cases  shown  in  Fig.  3  nd 

(d). 
With  the  simplified  rules  of  making  a  standard  block  di  .  .,: 
from  the  ird  determinant,  a  valuable  tool  is  developed  for  many 

purposes.    It  can  provide  a  student  with  a  means  to  db    v.    ind  visualize 
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a  cor  piex  feedback  system,  such  as  an  airframe  from  a  group  o:    !  Eferen 
tial  equations.    For  the  research  groups,  it  makes    malog  computei 
progr^rnin      ,  simplified  matter  with  respect  to  amplifier  usage.    In 

this  latter  case,  pitfalls  do  exist.    The  airframe  is  no  simple  mattei 

(a) 

The  reader  is  referred  to  Howe  uJ  for  an  excellent  coverage  of  the  air- 
frame problem  for  analog  computers.    Most  important,  though,  is  th< 
fact  that  this  standard  form  of  Chu°s  is  valid  for  any  multi-loop  lin<. 
feedback  system  whether  it  is  electrical ,  mechanical  or  a  combination 
of  the  two.    This  then  puts  a  regular  form  to  the  multi-loop  control 
system  from  which  design  or  synthesis  can  proceed  in  an  orderly  manner 
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The  Control  Problem  for  an  Airframe 


The  previous  three  sections  have  discussed  only  the  airfi 
itself.    It  is  quite  reasonable  to  say  that  while  these  point  out  m  my 
of  the  responses,  they  achieve  little  toward  solving  the  control  pro: 
In  simplified  terms,  the  control  problem  is  resolved  into  the  block   . 
gram  of  Fig,  4-1  „ 


i 


1 


-*■ 


Controller 


Controlled  System 


Figux 


A  S!  CONTROL  SY! 


Indeed  this  diagram  is  an  oversimplification  of  the  aircraft 
control  problem.    However,  it  does  point  up  the  fact  that  the  previc 
three  chapters  have  dealt  only  with  the  controlled  system  block.    Any 
complete  study  of  the  airframe  control  problem  must  include  the  otto 
block  components  of  the  amplifier  and  the  error  detector.    The  problei 
is  a  multiloop  type  and  of  no  small  magnitude. 

For  the  control  problem,  the  diagramming  becomes  a  simple 
extension  of  the  block  diagram  developed  in  Section  III,    Consider  the 
standard  determinant  for  the  longitudinal  system. 
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The  determinantal   an  iy  for  the  longitudinal  system  is  as  iollowo 
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From  this  a  block  diagram  can  be  made  up  from  the  determin  .r 
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This  diagram  is  of  the  three  degrees  of  freedom  longitudinal 
case.     Combining  this  to  a  control  system,  the  following  diagranr 
results  . 


45 


w 

CO 

o 

CO 


i 

S 

o 

i 


I 

0) 


# 

fo 


is,  the  control  problem  ior  the  simple  attitude  &     is 
presented  in  Fig.  4-2  and  3.    The  two  represent  the  different  point,  of 
view  of  the  aeronautical  and  the  electrical  engineer  although  the 
is  the  same.    In  Fig.  4-2  the  aeronautical  engineer  is  faced  with  t]  • 
primary  problem  of  compensation  within  the  airframe.    For  the  electrical 
engineer  who  inherits  the  fixed  airframe ,  the  attention  would  shift  to 
Fig.  4-3.    In  this  figure,  Gj  would  represent  the  pilot  response  in 
the  control  lever.    G2  would  represent  the  servo  system  response  o; 
hydraulic  or  electric  system.    Gg  would  represent  the  error  detection 
and  correction  signal  gained  from  an  electronic  element  of  some  \   n  . 
Similar  situations  of  block  diagramming  would  exist  for  the  changes  in 
speed  and  altitude.    For  the  most  part,  the  compensation  schemes  to 
achieve  stability  or  improve  stability  are  fixed  with  respect  to  j.  it! 
The  change  in  aerodynamic  coefficients  will  vary  the  paths  in  the  air- 
frame response  group.    If  the  electronic  devices  are  to  be  used,  then 
the  pick  offs  must  be  at  the  system  outputs  and  the  signals  fed  bad 
to  near  or  at  the  input  or  reference  signal  node  point. 

In  Fig.  4-3  the  pilot  constitutes  the  outermost  loop.    How*.:  ei 
this  is  not  n  mdatory.    In  a  fully  automatic  system,  the  pilot,  circu 
is  opened.    Few  paths  exist  where  both  aerodynamic  or  electronic 
means  could  be  used  to  achieve  stable  control .    One  can  be  noted  on 
Fig.  4-2.    That  is  elevator  position  response  71  which  could  be  alters  i 
aerodyn  .'   .c  illy  by  changing  parts  of  the  present  feedback  block  or 
adding  another  electronic  control  block,    It  is  then  fed  into  the  Ion  <_ 
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input  node  marked  Z^  Fe .    There  is  no  guarantee  that,  this  path  ^o.    et 
any  stability  problem  albeit  the  flutter  problem  may  be  aided.     It 
serves  as  an  example  of  two  possible  types  of  compensation  ovei  n  <_ 
same  path. 

The  lateral  systems  are  not  shown  or  discussed  in  this  chapter 
but  the  same  principles  of  control  systems  apply.    It  should  be  remem- 
bered that  the  airframe  response  shown  here  covers  the  small  pertur- 
bation response  concept  and  thus  is  greatly  simplified,,    The  following 
sections  will  deal  with  expansion  of  the  concept  for  some  workable 
systems . 
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Expansion  of  the  Standard  Form  for  Greater  Disturb  an 

The  system  developed  by  Chu  was  for  linear  functions.    It  is 

the  purpose  of  this  section  to  explore  the  equations  of  motion   ind 

determine  whether  the  equations  of  motion  can  be  used  for  largei 

turbances  than  considered  in  Section  III.    If  so,  then  to  what  de  jree 

can  this  system  be  used. 

For  all  practical  purposes,  the  small  angle  approxir     tion 

quite  well  up  to  one-tenth  radian .    Therefore  ,  8  ~      sin   (3  and 

cos^^r:  1  o    With  this  simplification  in  mind,  Equation  11-19  is  then 

rewritten . 

m  (U         tf  -  R7)   *  mz  sineQ  -  Wmg   (cose0  sin^Q)   + 

e  m    (cosQ„  cos  &  )   «  XA  +  X'U  +  X»W  +  X«W  +  X»n   +  X"  V 

+  XV)?    +  X*V  +  T     cos©      +  T     cos  A  mil  +  T  r  £  ,OSA 

tJ  I         V  o  T         u  "  Tu        x£  riw  £  RFM  coseT 


?=1 


The  thrust  forces  are  combined  into  the  forward  velocity  motion 

u  and  a  steady  state  reference  condition  of 

m  (n  W     -  R  V  )   +  m^  sine,,  s  X     +  T     cos©  „  V 

o  o         o  o  o         o  i 

is  assumed . 

Subtracting  V-2  from  V-l,  the  following  equation  results 

m  (u  +  Q„,„  -  +  ow  ■=■  R       -  V  r  «  vr) 

ov         o 


(cos©    sin<z5  )   +  &mz   (cos£^0  cosjZ$    ) 


4  j 


The  product  tern  qw  is  dropped  as  a  small  value  product  situ  ition 
The  sane  applies  to  the  vr  term. 

Converting  to  La  Place  notation  and  dividing  by  the  r;  iss  the 
following  equation  will  result 
si     *  0  ..'\  AX    +  sW    £}-'4/  IT   /6  -    ty '?.  cos  A    sin  (6 

+  6  *  ros  6     cos  JZ^    -  y  U  -  (sX*  *  yTT)  A  «*  -  sXq  © 

c    '     /   o  w         w  V-4 

[sT-J    +  x>^)7(  -  *#&      *  ° 

The  y  and  z  direction  force  equations  will  form  up  in  the  sane  in  mnei 
s  V  +  ^  oU  +  Uo  ^  "  ^  oTJo  **    -  WQ  0  -     <f  g  sin  aQ 

-  g   (cosAo  sin^0)   -  g^coseo  coS/0o    *  !#/? 

V-5 

-  sYp^  -  sYr    f-   (If  )£    -  (sY  £  +  Y£ )%     =0 

sU0A^     *  P0U0d  *  Vog  ^  -  Qou  -  "os  a  -  g  sinA0  6 
cos6>c  s^r£Zf0  <p    -  g  cos©     nosC6Q  »  Z^U 

-  (s  •"  0   U0  A<*  -  sZ    6  -   (s7t[    +  7  rj)7|    -  0 


'.0 


The  moment  equations  can  be  formed  in  much  the  same  manner.     From 
Section  Two,  Equation  11-22  would  convert  to  this  equation: 

sy  ■     </.  »  SH0a)  -  e    (*p0e-  sqo^)  (^ 


In  a  similar  manner  Equations  11-23  and  11-24  can  be  brought  to  a 
similar  maximum  linear  condition.    They  are  shown  in  the  following 


equations: 


*(t£.Ys'Po&  +  sQo^    +  -L    (j5?oft  "   Rpo<//) 


,,)  !J    ^o(«s   (s,r-   -  F  )©       (sM^    *      »  ^X, 

V-8 


s2^  f/— l(sPoe  +  B°o^    f  "     (sQo  ^+  sPoe  "  SV) 

/S  ^S~  sNR</>  -  Nr|      -   (sN^  4   N£)   ^  v-9 


- 

is 


Also  the  control  hinge  moment  equations  III-9,  III- 31  and 
111-32  are  recalled  in  their  same  form. 
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These  nine  equations  represent  the  linear  equations  of  motion 
of  the  airframe  and  include  all  cross-coupling  effects.    These  equations 
still  retain  the  requirement  of  a  linear  approximation.    The  equations 
with  the  author's  modification  of  adding  the  equations  -1  +  1  =  0  m   -  e 
up  Table  7-1.    This  is  an  18  x  18  determinant  .    From  this  table  the 
block  diagram  of  rig,  5-1  can  be  drawn . 
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TABLE  5-1    DETERMINANT  ARRAY  OF  AIRFRAME  WITH  LINEAR  APPROXIMATIONS  FOR  SEC  DEGREES    ^  FREED0M 


FIGURE  5-1     BLOCK  DIAGRAM  OF  THE  AIRFRAME  AND  CONTROL  SYSTEM  USING 
THE  STANDARD  FORM 
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block  diagram  then  represents  the  complete  airfi      e 
respon.  e  m  lintaining  linear  approximations.    Even  then  thj  ei 

t  ons.    As  an  example,  consider  some  airplane  flying  with 
a  steady  state  wing  angle  of  attack  of  15  degrees,     (assume  a  "con- 
vention" wing  and  the  reference  axis  is  the  zero  lift  line)     Any  pilot 
will  quickly  testify  that  an  increase  of  5  degrees  angle  of  attach    tl 
that  point  will  probably  result  in  a  rather  wildly  oscillatin  j  -    meu   er 
popularly  known  as  a  spin.    This  condition  is  generally  outside  of  tl  e 
range  of  the  aerodynamic  forces  and  moments  usually  computed  „    So 
aerodynamic  coefficients  such  as  lift  (Cl)  are  decidedly  non-linear  in 
this  region.    Therefore,  the  restrictions  imposed  on  this  derivation  of 
terms  must  include  a  restriction  that  the  airframe  is  not  operating     I 
the  very  edge  of  its  flight  envelope. 

The  system  has  been  diagrammed  directly  from  Table  Y~l„    1   e 
terms  in  Column  One  are  all  connected  at  node  I  of  Fig„  5-1.    Likev. . 
the  terms  in  Column  Two  of  Table  V-I  are  all  connected  to  no<.e  11  ,    E> 
nodes  appear  in  the  beginning  and  end  of  the  block  diagram,    rhe   e  ire 
inserted  as  provisions  for  control  functions.    These  control  nodes  are 
lettered  alphabetically  .    The  block  diagram  contains  the  required  loop 
for  aileron    md  elevator  control.    The  rudder  control  loop  can  be  e  isily 
inserted  when  desired. 

lx  iminatlon  of  Fig.  5-1  reveals  many  interesting  points  here- 
tofore known  but  sometimes  only  by  intuitive  reasoning  e  The  concept 
oi     ej.  ii  ite  :>y  terns  for  the  symmetric  axis  and  asymmetric  ixis 
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3  ;s  very  well  shown.    The  level  flight  steady  state  conditions 
used  m  Section  Three  stated  that 

Po  =  Qo  =  Ro  =  Vo  =  Wo  =  0 
ning  the  paths  between  nodes  VIII  and  IX  re v«.    ,  tl 
o-  the  perform  ince  functions  are  zero  under  these  specified  con  d . 
of  steady  state  flight,.    There  is  no  direct,  feedforward  or  fee  .      cl 
paths  between  nodes  VIII  and  IX.    Thus  the  restricted  case  of  steady 
51  tte  level  flight  conditions  with  small  perturbation  is  clearly  di    I   € 
into  the  two  systems      The  longitudinal  system  includes  the  forvv 
velocity,  vertical  velocity,  pitching  motion  and  elevator  position. 
The  block  diagram  for  the  system  can  be  further  restricted  to  two 
de  -rees  of  freedom  by  holding  the  X  velocity  change,  u,  at  zero, 
This  allow     (  h  mge  in  vertical  motion  and  pitch,,    This  would  then  b« 
Figs  -nd  3-7  of  Section  III. 

The  so-called  lateral  airframe  response  system  comprises  tl 
nodes  trc    IX  to  XVIII.    The  rudder  and  aileron  motion,    sideslip,       i\\ 
and  roll  comprise  the  motion  for  this  sytem.    With  the  steady  state 
level  flight  conditions  previously  mentioned,  the  block  dia  p    r-   is 
identic    1  : .,  }     .     \-h  of  Section  III. 

An  interesting  note  of  these  separate  systems  is  the  cros^ 
coupling  terma  that  begin  to  arise  with  larger  disturbances.    An 
portant  one  is  the  pitch  due  to  sideslip.    The  effects  are  labeled  M 
in  the  diagram.    Several  aircraft  are  subject  to  this.    Usually,  tl 
is  -ng  or  nose  down  motion  that  appears  as     ■  .      Lip  an  fie       8    ) 
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inert,  i   es,     in  the  simplified  systei  s  this  effect  wis  neglecte  i 
is  of  ne  ilue  for  very  small  disturbances,  but  can  be< 

notic  oy  the  time  the  sideslip  angle  (/•$    )  re^c she     £iv< 

.  te  effect  exists  and  is  shown  by  this  expanded  system. 

excellent  methods  of  analysis  are  available  for  these    ;1    n 

forms  brought  forth  in  this  work.    Methods  of  analyzing  outpu:      n 

(fn  (ill 
response  to  one  or  more  inputs  are  available.    Work  by  Thaler     """  eon- 

l    ins  a  very  useable  system. 

Where       fV      S     the  output  signal  at  the  n  th  node 

^      s    characteristic  determinant 

^    the  cofactor  of  the  determinant  with  BCa" 
the  node  location  receiving  the  inpul 
and  30nto  the  node  location  of  the  o::t; 
signal . 

A 

h\  —     the  performance  function  on  the  lnp    I 

prior  to  entering  the  input  node . 


1      s.     the  input  signal 


G^y  the  direct  path  performance  function  between 

the  output  node  and  the  output  signal  „ 

Thus  a  single  input  signal  would  yield  an  output  signal 

11 


A 

I 


For  two  or  more  input  signals  additional  subscripts  would  i  e  rnd 

the  functions  are  additive.    The  rules  of  superposition  of  linear  o  itpul 
justify  this  method. 

a  ■:.  ■     1 
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i  would  give  the  resulting  output  of  one  or      ore    nputa 
at  various  nodes  in  the  system. 

The  iegrees  of  freedom  can  be  restricted  also.    In  the  block. 
.     consider  a  cut  rrmde  in  the  path  at  point  "S". 


DIAGR 
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opening  in  the  path  now  eliminates  the        rial  to  the  d    t  < 
p    :       the  feedback  and  the  feedforward  path,    The  result  o:  ; 
determinant  Is  the  insertion  of  zeros  in  all  the  performance  function.? 
in  the  column.    It  leaves  only  the  constant,  one,  in  the  i       n  onal 

block,    The  reader  will  recall  that  the  main  diagonal  has  the  te 
1  +  PF     '       o    This  is  shown  below  as  opening  the  output  path  from 
third  node. 

x       x        0       x       x       x 

x       x        0        x       x       x 

XX  XXX 

X  X  0  X  X  X 

X  X  0  X  X  X 

X  X  0  X  X  X 

The  output  p-ith  from  node  three  is  opened  at  point  s  and  the  deter- 
minant Is  changed  in  this  manner. 

A  special  case  occurs  for  the  diagram  illustrated  in  I 
The  second  node  has  only  one  input  signal.    Because  the  path  v.  is  c  u.1 
at  point  "s'\  the  output  from  the  second  node  will  also  be  z&ro ...    Tl 
all  feedforward  values  from  this  node  will  be  zero.    The  direct,  an  . 
feedback  paths  of  the  second  node  are  non-existant  to  begin  \ 
rerrain  so.    The  determinant  form  would  appear  as  now  shewn 


x       x        0        0        x       x 

x       x        0        0       x       x 

x       x        1        0        x       x 

x       x        0        1        x       x 

x       x        0        0        x       x 

x       x        0        0        x       x 
Simple  manipulation  now  reduces  the  six  by  six  determinant  to  a 
four  by  four  determinant,    Therefore,  when  it  is  desired  to  consider 
a  problem  with  one  or  more  degrees  of  freedom  held  at  zero,  the  simple 
method  of  opening  the  output  signal  path  from  the  node  and  placing  a 
zero  for  all  the  performance  functions  in  the  column  leaving  only  the 
constant  value  one  in  the  main  diagonal  point  accomplishes  this 

Referring  to  Section  II ,  the  reader  can  observe  the  general 
equations  of  motion  for  mass  particles  and  examine  one  such  as 
Equation  11-22  which  is  repeated  below. 


'    '     ^K],    '  '  3-B)   -  E  (PQ  +  R)   +  c  K,Q  -  bK2R  - 

» 
L„  Uv  J   L'p  +  L'r  +  IA  i       +  JJ   t  V- 


%\   -  L't?    "° 


Note  the  third  term  QR  is  a  produce  of  (Q    +    q)  (R       «   r  ) .    The  pro- 

o  o 

duct  when  multiplied  is  Q  R         Q  r  •■-  R       +  qr.    From  the  equations 
of  Section  I: 

Q  *  6  cos  <jf>  +  <fsin  0cos  ©  V~ 
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(p      j  pc   s©-   ©sin   $ 
When  the  small  angle  approximation  is  no  longer  valid,  the  hand] 
of  this  product  term  becomes  rather  involved,,    Furthermore,  the  aero- 
dynamic coefficients  of  the  L^,  L '  ,  L'  etc  terms  lose  their  linear 
approximations  over  large  changes  in  flight  conditions.    Most  aero- 
dynamic coefficients  are  of  a  decidedly  non-linear  nature  over  the 
full  span  of  an  aircraft's  maneuvering  capability.    The  coefficients 
are  definitely  non-linear  with  Mach  number  when  the  velocity  range 
covers  the  subsonic,  transonic  and  supersonic  spectrum.    The  reader 
can  observe  the  general  equations  of  motion  in  Section  II.    Considering 
the  complexities  of  incorporating  equations  V-14  and  V-15  for  large 
movements  and  the  non  linearities  of  the  aerodynamic  coefficients, 
it  is  not  long  before  one  realizes  that  the  non  linearities  are  numero 
and  cover  all  the  equations.    The  prospect  of  block  diagramming  for 
the  purpose  of  placing  all  non  linearities  in  one  block  to  be  handled  by 
a  describing  function  scheme  is  hopeless. 

The  possibility  of  doing  the  problem  with  an  analog  computer 
and  using  no  approximations  is  fair.    Undoubtedly  the  task,  if  attempted 
will  be  of  considerable  magnitude. 

]  bus,  it  is  seen  that  this  system  of  standard  determinant  arrays 
and  standard  block  diagrams  can  be  used  to  handle  the  airframe  problc 
as  long  as  a  linear  approximation  is  maintained.    The  system  is  flexible 
and  can  c  irry  varied  conditions  of  multiple  inputs.    It  can  also  be  note 
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that  for  the  six  degrees  of  freedom  problem  over  a  large  variation 
it  is  of  limited  value  except  for  diagramming  a  computer  problem. 
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6o  j  rospect  and  Retrospect 

The  airframe  motion  problem  has  been  derived  quite  genei 
in  the  first  two  sections  of  this  work.    The  attempt  was  made  to  mainl 
as  many  effe<  ts    is  po  s    ;ble  in  the  general  equations  of  motion.    This 
was  done  so  that  the  basic  equations  of  motion  could  be  used  in  a 
variety  of  airframes  varying  from  missiles  to  ground  effect        <  'lines. 
In  the  succeeding   sections  the  author  has  confined  his  discussions  to 
equations  of  basic  airplane  types.    Rotor  effects  were  disregarded  and 
thrust  effects  were  combined  with  drag  to  determine  the  effects  along 
the  X  axis  system ,    This  should  not  imply  that  the  equations  of  motion 
cannot  be  used  for  other  airframe  types.    They  can  in  all  cases  where 
the  assumptions  of  a  rigid  body,  constant  mass,  etc.  are  still  valid. 

The  author  has  several  thoughts  on  areas  of  further  exploration. 
He  will  discuss  them  briefly  and  mention  any  work  of  possible  interest 
in  connection  with  this  field. 

The  primary  aim  for  control  engineering  in  this  field  is  the 
development  of  an  orderly  process  for  multiloop  compensation  „    It 
shown  clearly  in  Fig„  5-1  that  several  paths  of  possible  compensation 
exist.    What  path  or  paths  are  the  most  successful  and/or  easiest  to 
achieve  multiloop  compensation?    Some  work  has  been  done  in  this 

field.    Specifically,  an  electronic  compensation  of  a  missile  using  t> 

i  1 
determinant  method  was  done  by  Anderson  and  Roane.        They  confined 

their  efforts  to  the  roll  and  yaw  coupling  problem  for  a  specific  vehicle . 

Their  efforts  achieved  a  satisfactory  solution  to  meeting  the      I      ile 


specifications,  but  no  definite  method  of  analysis  as  to  the  best  or 
poorest  i-   ";   to  compensate  a  multiloop  system  was  generated.    A 
possible  approach  would  be  an  examination  of  aerodynamic  coefficients 
and  their  known  effects  .    In  chapter  ten  of  Perkins  and  Huge    there  is 
a  discussion  of  the  fact  that  drag  or  coefficient  of  dray  markedly  affec  I  . 
the  damping  in  the  longitudinal  phugoid  motion  of  an  airplane .    The 
airframe  short  period  is  affected  largely  by  the  elevator  hinge  moment 
coefficient-.  C  and  C  .    An  examination  of  these  factors  an  ! 

others  might  lead  to  a  pattern  which  could  successfully  solve  the 
problem  for  any  multiloop  system. 

The  block  diagram  development  should  provide  a  basis  for 
analog  simulation  of  the  linear  problem  expressed  in  sections  III  and 
V.    The  analog  computer  arrangement  may  be  further  expanded  to  incl 
some  of  the  non-linear  effects  using  function  generators,    This  proLle: 
allowinj  six  degrees  of  freedom,  would  overtax  the  capacity  of  tl 
analog  equipment  presently  available  at  this  institution  „    If  an     r.   log 
computer  with  a  100  amplifier  capacity  and  50  function  generators  was 
provided,  the  problem  could  be  nicely  handled 0    Most  assuredly,  co 
binations  can  be  made  by  grouping  various  functions  which  will  cut  the 
above  numbers  by  one  half.    However,  various  outputs  will  not  be 
available  for  analysis  and  many  signals  which  are  not  physically 
realizable  in  the  actual  system  will  appear  as  outputs . 

The  high  speed  digital  computers  could  be  utilized  in  a  project 
of  statistical  evaluation  of  compensation  schemes.    The  co 


•  4 


c  ip    \    Lity  of  producing  hundreds  of  solutions  per  minute  would  allow 
ex   -  In  ttion  oi  data  to  determine  trends  and  patterns  for  i  or 

pen: "-.tor  schemes.    A  successful  compensation  criteria  i        t  be 
derived  in  this    n    nner. 

An  aerodynamics  group    ■     /  examine  four  items  of  interest  in 
this  field.    The  control  surface  hinge  moment  equation     use     in  this 
work  were  extremely  simplified  „    Exploration  of  these  effects  to 
attain  more  refinement  would  be  justified  „ 

In  section  II ,  the  author  pointed  out  the  absence  of  compu- 
tational methods  for  various  stability  derivatives  „    Further,  the 
assumption  of  a  quasi-steady  state  condition  is  made  in  corap  i1  In  • 
the  known  derivatives.    A  study  by  Etkin     ^  cites  an  improved  S}  >te 
ro  attain  better  accuracy  of  the  aerodynamic  stability  deriv  ite 
Further  work  in  this  area  appears  justified . 

There  is  also  a  need  to  determine  the  effects  of  aeroe!      tii    I 
on  the  problem.    To  this  end,  instrumentation  and  flight  testing  an 
airframe  might  prove  useful.    The  capability  of  the  Naval  Air  Facility 
to  r    mi  .in    md  oper  ite  many  types  of  military  aircraft  should  permit 
use  o:     or  e  airplane  or  aircraft  for  this  work.    An  additional  use  wo 
checking  the  effects  of  automatic  stabilization  equipment  in  use  on  tl  e 
airframe.    The  HSS-i  helicopter  has  an  excellent  control  system  to. 
this  type  of  research.    However,  the  equations  of  motion  used  in  thj 
work  will  require  some  refinement  to  control  the  \  .  inertia  problem 

c  iu  3©  I  by  the  tilting  rotor  disk. 
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A  final  topic  for  consideration  is  that  of  signal  flow  theory, 
This  idea,  which  has  been  advanced  in  recent  years,  may  be  com- 
patible to  the  same  "standardized"  system  as  presented  here  for  ti 
block  diagram  method . 
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7.  Conclusion 

is  been  conclusively  shown  that  Ghu's  ,Dst  ,n  1  rd"  bloc1 
diagram  and  determinant  method  can  be  applied  to  a  linear  system 
other  than  a  servo  control  system.    With  some  modification 0  any  set 
of  linear  differential  equations  can  be  block  diagrammed  into  a  stand 
form  for  general  analysis,    A  very  workable  system  exists  for  analysis 
of  a  multi-loop  control  dystem  using  one  or  more  input  signals,,    The 
airframe  stability  problem  can  be  examined  for  the  linear  approximation 
using  this  very  effective  system. 


67 


REFERENCES  AND  BIBLIOGRAPHY 

1 .  Chi: ,  Yaohan;  ""A  Generalized  Theory  of  Linear  Multi-loop 
Auton  itic  Control  Systems",  Fh„D  Thesis,  MIT,  May 

2.  Dooiin,  Bo  F„;  "The  Application  of  Matrix  Methods  to  Coord  n 
Transforms  Occuring  in  Systems  Studies  Involving  Large  Motions 
on  Aircraft"  ,  NACA  Technical  Note  3968  ,  May 

3.  "Dynamics  of  the  Airframe",  BuAer  Report  AE-61-4II  prepared 

the  Servomechanisms  Section  and  Aerodynamics  Section,  Northrop 
Aircraft,  Inc.,  195  J  . 

4.  Etkin,  Bernard;  Dynamics  of  Flight,  Stability  and  Control, 
John  Wiley  and  Sons  Inc.  ,  New  York,  1958. 

5.  Perkins,  CD,  and  Hage,  R„  E.;  Airplane  Performance,  Stability, 
and  Control ,  John  Wiley  and  Sons  Inc.  ,  New  York,  1349  „.. 

6.  Duncan,  Wc  J.;  The  Principles  of  the  Control  and  Stability  of 
Aircraft ,  Cambridge  at  the  University  Press,  1952c 

7.  Charters,  A,  C;  "The  Linearized  Equations  of  Motion  Under- 
lying the  Dynamic  Stability  of  Aircraft,  Spinning  Projectiles, 

ni  Symmetrical  Missiles",  NACA  Technical  Note  3350, 
January  1955 . 

8.  Howe,  Ro  M.;  "Coordinate  Systems  for  Solving  the  Three 

Dimensional  Flight  Equations",  WADC  TN  55=747,  June  ] 

9.  "Dynamics  of  the  Airframe"  op  cit 

10.  Thaler,  G»  J.;  "Linear  and  Non 'Linear  Compensation  Theory", 
Prepared  Notes  for  Course  EE  676,  U.  S.  Naval  Postgraduate 
School,  Monterey,  California  1961,  p  19-21. 

11.  Anderson,  R0  W    and  Roane,  D.  P.;  "An  Investigation  of  the 
Roll-Yaw  Couplings  in  a  Surface  -To-Air  Guided  Missile", 
Masters  Thesis,  MIT,  May  1960. 

12.  Etkin,  Bernard;  "Aerodynamic  Transfer  Functions:    An  Improve- 
rent  On  Stability  Derivatives  for  Unsteady  Flight00,  UTIA 
Report  No  42,  Toronto,  Canada,  October  1956. 


■  b 


13,  Smith,   1       •na  Hick.,  W.D.T.;  "The  R.A.E.  Electronic  S  itor 
for  Fl   ttei  Investigations  in  Six  Degrees  of  Freedom  or  L« 
Aeronautical  Research  Council,  Great  Britain ,  R„     n  .  M  ,. 
Sept..  •  bei 

14.  :  irov,  7.  Au  and  others;  "The  Airplane  as  an  Object  of  Control", 

(Translation  from  Russian) ,  NASA  Technical  Transl  ition  I 

Was:    n  .'on,  D,C,  ,  October  1959. 

15  o  Winchester,  F.  W.;  Aerodonetics ,  Do  Van  Nostrand  Company, 

New  York  1909 . 

18.  Durand,  W.  F„  and  Jones,  B„M.,°  Aerodynamic  Theory,  Volume 
"Dynamics  of  the  Airframe",  Durand  Reprinting  Committee,, 

California  Institute  of  Technology,   1943. 

17.  Abzug,  M.  J.;  ""Kinematics  and  Dynamics  of  Fully-Maneuverii 

Airplanes",  Douglas  Aircraft  Company  (El  Segundo) ,  Santa  Monica 
( .  alifornia  ,  June  1 9  5  2 . 

18.  Wykes,  J.  H,  and  others.  "An  Analytical  Study  of  the  Bynami< 
of  Spinning  Aircraft",  Part  One,  North  Ameixan  -V  -ion  Inc=  , 
December  1958. 


I 


Al  .  ..ND1X  1 

BASIC  EQUATIONS  OF  MOTION 
i    .     ;ic  c  oncept    of  mass  in  motion  and  the  equations  of 
■  Ion  c  in  he  derived  in  -vector  notation.    This  is  done  using 
particles  in  motion  and  determining  their  effects  in  two  separate  axj 
systems.    The  basic  reference  system  shown  in  Fig.  AI-I  is  th« 
axis.    This  is  the  non-rotating  reference  system  and  is  also  known    •. s 
the  Galilean  Axis.    The  axi      -ystem  is  fixed  in  the  rigid  mass  but  ; 
rotate  about  the  fixed  space  axis  system.    It  is  the  relative  system. 
A  special  set  of  notation  is  used  here  in  this  appendix  which  differs 
from  the  body  of  this  work. 

The  definition  of  the  special  terms  in  this  appendix  is 
follows* 

i  =       the  vectorial  unit  ilong  the  X  axis 

—>  4 

j  =       the  vectorial  unit  along  the  Y  axis 

-*• 
k  =       the  vectorial  unit  along  the  Z  axis 

#  &  * 

i  '  =       the  vectorial  unit  along  the  X    axis 

******    & 

j  -        the  vectorial  unit  along  the  Y'     axis 

k  the  vectorial  unit  along  the  Z    axis 

0  ~        the  on-jm  of  the  relative  axis  system 

*  4 

0  =       the  origin  of  the  space  axis  system 

R  =       the  vectorial  distance  from  the  space  axis  origin  to  the 

relative  axis  origin 

r  =       the  vectorial  distance  from  the  space  axis  origin  to   i 

i  tide  in  space 
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P  =       the  vectorial  distance  from  the  relative  axis  origin  to 

the  mass  particle 

A 

dm  =       the  mass  particle 


t 


=       the  vectorial  distance  from  the  mass  center  to  the  origin 
of  the  relative  axis  system 


J\  =       rotational  velocity  vector  of  the  relative  axis  system  abc 

the  space  axis;  an  absolute  velocity  vector 

-*  A 

<jO  -       absolute  rotational  velocity  vector  of  a  rotor  whose  axis  is 

attached  to  XYZ  axis  system 


I: 


momentum  vector  for  the  entire  mass  body 


H  =       moment  of  momentum  vector  with  respect  to  the  space  axis 


-> 


A 


G  =       moment  vector  with  respect  to  the  relative  axis 


-> 


A 


F  =       the  summation  of  the  external  forces  on  a  body 

/  ft       the  summation  of  the  interna!  forces  on  a  body 

This  appendix  will  derive  the  basic  equations  of  motion  for  any 
mass  moving  in  space.    After  the  basic  equations  have  been  developed, 
the  effects  of  control  actions  will  be  added.    The  final  equations  will 
then  be  converted  to  standard  aeronautical  terminology  which  is  used  in 
the  body  of  the  work  for  aircraft  stability  considerations , 

These  equations  have  no  doubt  been  derived  several  times  in  the 
past  one  hundred  years.    However,  the  author  was  unable  to  find  any 
clear  derivations  in  any  of  the  reference  material  covered,,    Therefore 
the  derivations  in  this  Appendix  are  submitted  for  the  reader's  interest 
and  for  background  information . 
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ody 


Consider 


->       -*       -» 

f*  s  ix  +  ,1y  +  kz 


Momenturn  vector  for  the  whole  body  L 


=  J  dm  dr_ 

J  dt 


dr_  =    (Vo   ♦   (u)x/°)    *  f* of  roi 

dt 


or) 


L  =  J  -  /dm   (Vo  +   ftO"^    +  /°  rotor) 


dt 


(  r      -f 

J  J 


■  ++  F 


but  summation  of  internal  forces  equals  zero 


Hence 


dL     » 

dt  y 


d   r 


2*      ~> 
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For  an  airframe  with  a  large  rotor  (ie  helicopter ,  STOL  aircraft  or  GEM) 

consider  addition  of  the  rotor  to  the  airframe  using  the  mass  center  oi 

the  aircraft  as  the  origin  of  the  axis  system. 
c\?T         d2R  +A_x  Z'+Jlx   (ax   />),  +  /*+  2/Lx  /> 


dt?       U- 


J 


i     isport   terr.s  relative  coriol:' 

acceleration     effect 
term 


Co  ■■. 

;-•  -  rt  y  f+£L*  (Ax  /°)  -  />+  ?d  r./° 


o 


I 

•     ■  • 


."'TT-    -   i-;:-',v  +   K"-;'-v?.  +  i-"-( 


(rt-.Vz   --fl-JTy 


)  +  .1*  (il3    -  -il^ 


k^UL^y  -  il2Vx)    +0  -  0  +  0+  0  whjjfre/*  -  0 


m  d 


at* 


v\ 


v-     il/~  -XV, -^ 


=  j*(vy  H-n3vv    A*,V,)m 


ul  3' 


k*Fs   =  k*(Vz   +AjVy  -A2Vx)m 


^ 


neglecting  rotor  effect. 


Therefore  in  scalar  forw 


: 


7     +A?  v   -A3Vy)m 
(Vz   4-A^y  -A2Vx)i 


Concurs  with  Etkin 

lwU.3 

/lso  Duncan  Page  6U 


Figure  AI-2     VECTOR  DIAGRAM  OF 
ROTOR  HOT ION 


PQ   is  fixed  in  the  rigid  body9     The  rotor  spins  at  an  abso. 
angular  velocity  of     but  the  mass  center  of  the  rotor   system  is  at 
point  Po     Therefore ,  ftc  and^c     =  0e     It  is  also  assumed  that 
_/TL  =  constant*     Therefore,  the  additional  term_fLx(.n,x/ac) 
enters  the  problem, 

i        j        k 


-0.   X    (°:    ■ 


x    '    y       a 


A  x  (3  x  G;   =   (A» 0)3  -  (A*B)C  is  a  basic  identity,  therefore 
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Note:      coordinates  i ,  j ,  k  are  those  in  the  body  and  not  the  inertia! 
or  c-  irtl    system.    The  values  of  -A,  are  primed  for  the  same 
reason.    It  is  quite  possible  that  /^       will  be  0  considering 
alignment  of  the  rotor  on  the  xz  plane . 


I  .   -  F,  +  F«  -  m„ a  +  m^a     F,  s  airframe       Tn  =  rotor 
r        1         £•         i  *  1  I 

In  resolving  this,  it  must  be  noted  that  the  mass  centers  differ 
from  those  of  the  combined  masses. 

It  should  be  noted  that  the  equations  for  the  acceleration  o! 

particle  also  contain  rotational  and  relative  velocity  terms,  however, 

/     and  /^  =  0  as  these  particles  are  in  a  rigid  body.    Further,  the 

summation  of  particles  about  the  mass  center  makes 
"f      -*       ->■  -*       -*■ 

c\m    [_( JL  X  P)  1r  J)_    X   (-^-  XPJJ 

Moment  terms  may  be  considered  using  Moment  of  Momenta; 

From  this  the  differential  may  be  taken  and  the  moment  is  then 

derived.    Consider  H     to  be  the  moment  of  momentum  of  the  system  with 

respect  to  the  inertial  (non-rotating  axis) .    The  derivation  will  include 

a  rotor  spinning  on  the  rigid  body . 


Figure  AI-3  MASS  DIAGRAM 
FOR  MOMENT  ANALYSIS 


M  =  mass  of  the  rigid 
body  less  rotor 

m  =  mass  of  the  rot 


The  absolute  angular  velocity  of  the  axis  OXYZ 


is  fL  =  •  /V,  +  j  ^i-^Aj 


->     ~Z 


-> 


The  absolute  angular  velocity  of  the  body  is  uJ  -  i  u;,  -*-j  ^+  K  ^5" 

The  relative  angular  velocity  of  the  body  is  J ^   -Yi-/  =  n  feco 

The  angular  motion  of  a  rotor  in  motion ,  but  attached  to  the  body  will  be 


UJ 


-jtWO,  «/Ll^7Cu:z-AL5+f?(u;:}-/l.3'). 


-* 


Henceforth  this  relative  term  will  be  noted  as  K  (KAPPA). 

* 
The  equation  for  H     is: 

-  /r  x  dr  dm  where  r  3  P.  +  /* 


Recalling  that  the  momentum  equation  was  L  = 


-> 
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then 
->       -*      -9 


■j* 


-> 


t 


--?* 


The  term  under  the  integral  will  be  known  as  H 


Recalling  that      ^    ■  Vc   +/lx /**  +  /*<      t      ) 


-^  -^     ■-* 


H=//?        rv-»       'Ax/"    '  /°(ro1    j 


This  can  further  be  resolved  into  the  moment  of  momentum  about  its 
own  Oxyz  axis  and  the  velocity  moment  with  respect  to  the  inertial 


axis  . 

F  -  (M  +  m)/^  x  Vo  */(/*x/Lx/^  dm  +  //*  x   (l(  x/3)  d 
Hb  =//^x  (itx^  din 


m 


To  derive  moments  consider  recalling  H     =  /    (i 


-5» 


—-  x  £1 
rlt      at 


)*■: 


(F   +Z3)  x  c,?r 
dt? 


-1 


rid  rec?ll         tl    t  v  «  dl  -I  f  r     on 

dt     I       ^ 


Then  substituting  these  in, 

-»_►-»->'  -*  A 

T    +  G  w'.i 

dt 


A 


J 


However ,  we  also  noted  that 

:  ]     •  /  (/  -  dr   ]dT 
dt/ 
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Therefore: 

oH*     «  dE     -  L  +  R  -  dL         dH 
dt  dt  dt 


Equating  the  two  results  of     ^H*   it  is  seen  that 

dt" 
-*      ->      ~->    ->     v       -^       •> 

JNtfL     +  G  -  dR  x  L  +  R^S^L  +  £12 

d&  dt  eX     dt 

Recall  that     .     -     . 
-^     — ^ 

dR  *-  Vo 
dt 


and  that 

1  -  m)(Vo  l/bc/J.  +  /%  rotor) 

i  -  d       |(M  +  m)/~.  x  Vo  *  Ho  +  Fo^t 

dt  L  ^1J 


[H     «  (m  +  M)     d/9      3    *ro  +/*   x  dVc 
dt  U~ 


dt 


T'V    +  Horel 


combinj 

G  *  Vo  k   (Vo   'Arrp)  T'I  +  Vo  x  (7o  +Jlx/        +  />     rotor)  m 


— ^  —  ■  -*^ 

+  (M  ♦  m)  /°„  x  dVo     +  d       I  w0  +  !  q     ,  1 
dt  dt     L  J 


,  x  Vo  +  (ilx/5  )  x  Vo)   M 


A  further  assumption  is  stated  at  this  point.    This  is  that  the 
point  of  rotation  of  the  rotor  is  fixed  on  the  body  and  remains  so0    The 
shaft  does  not  have  nutation  about  some  axis  but  remains  fixed  on  a 
mount  of  the  rigid  body,    Therefore,  the  /";        of  the  rotating  body  is 
zero.    Also,  logically,  in  a  rigid  body  the  mass  does  not  shift  hence 


f« 


is  also  zero 
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In  crossing  a  vector  with  itself  Vo  x  Vo  the  result  is  zero. 


.  O    Y.     70    ~ 


i         j         k 


-r  Tr 

'1        "2         3 


-  a    (VpV3  -  V2V3)    <    i  (V^  -  VjVj) 


+  k   (V]V2  ~  Vf/3)   ^  0 


Vo  x  Tro  ■  0 


Furthermore ,  crossing  two  dissimilar  vectors  in  opposite  order  yields 
opposite  signs  and  equal  magnitude. 


70   Vo 


v 
2         3 


A   /°?    /^ 


-» 


i  (v2/>3-/y3)  *  j  (/^-/Vi* 


hs>    -^ 


/°c  x  Vo  -  i  (/*2v3  -  /»3v2)  4  7(^3-, -  /^V  )   *  k  ^V2  -  r 


-*•       — > 


Thus  Vo  x     Z5     +  /*    x  ^o  -  0 


Incorporating  these  two  facts,  the  equation  for  G  can  be  simplified  to 


the  following: 


|_  p*^relj+  <H  ♦■).</* 


x  dVo) 
dt 


— >  -^  ->     -* 

r)-.'  •  a  n0  -/]_x  ji0  -;n  a  general  enuation. 

dt 
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-7  -* 

Thus  it  is  necessary  to  solve  dHo    to  find  G„ 

dt 


-> 


Considering  Ho  =  j  [_/     /  (c<J   x  l)\  which  accounts  for  the 

absolute  '  elocity,, 

Recalling  an  identity  of  vector  algebra: 

.  x   (p  xc)  a  (A*C)B  -  (/.•B)cthen 

Ho  -y    (x2  *  y2  *  OCiL^  +  TtJ2   +fu)3)   -  (xc^  +  yu>?  4  z  cJ3) 
(ix  +  Jjr  *  k?)       dn 


<■  k  J^?  +  y;'u>3  -  xs  ^  -  yzu'2l 


dm 


The  above  terms  can  be  recognized  in  part  as  moment  and  product  of 


inertia  terms . 

Ix 

•)dm 

iy 

=*  f'J  «  z' 

Iz 

&fh?  *r 

')dm 

Ixy 

=  J (xy)  dm 

Iyz 

=  J  ' .  -  )d  t 

Ixz 

=  J: 
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Therefore: 
Ho    -   I    (Ix^,     -   jv-u: 


-'  ~  J 


oczcJ^:    •   j   (-IxyiA 


Tyt 


0U-, 


Jyzi 


+  k  (I*  u}7  -  X*z  iA|_  -  Iyz  u>2) 

Considering  the  effects  of  the  aircraft  as  a  rigid  body  first 

dHo  e  Ho  *ilx  Ho  c.nd  for  a  rifid  body  u)  -  -1  J- 
dt 

1  (TxTLj    -  L:y/{2  .  Ixzi.3)    +  J   (-Eg^  +  I7/L2    ■  V-^3) 


*  k  (-Ley 0. 


Iyzil2  +  Iz./13)  + 


i 

Al 


A 


JX. 


(  IV = Ty^  ~  T"  " )      ( Ty— Txv«=Tvz ) 

- 


For  any  ordinary  aircraft  or  missile,  the  product  of  inertia  terms  Ixy 
and  Iyz  can  be  considered  zero  because  of  the  mirror  symmetry  aspect, 
Ixz  may  be  zero  if  the  principal  axis  is  selected.    Thus  with  no  los:  In 
generality,  consider  Ixy  =  Ixz  =  Iyz  =  0 


dHo  -':     [pcOj    +  I7J12^3  -  Ty-/l?/lJ 

+  ?JVd2  -  Bcfl^  -  --^1^3] 


However ,  including  the  effects  of  Ixz,  the  product  of  inertia  term, 


which  is  often  not  zero 

s  - 
dt  ' 


4 


Using  standard  aerodynamic  terminology,  the  following  items  are 

defined: 

A  =  Ix 

B^Iy 

C  =  Iz 

E  =  Ixz 

Q  =A2 

RSn3    ' 

Thus 

G  -    dHo     where  system  is  chosen  such  that  /       *  o 

dt*  C 


i   ]a  /**-  ER  -  EFQ  +  QR  (C-B)J 


If." 


j    [BO   +  PR  (A-C)   +  E  (P?  -  R?)] 

k  fcR  «  tt*  i  r    (b-c)  +  eqrj 


+ 
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These  equations  concur  with  Duncan,  Page  66,  and  Etkin,  Chapter  4„ 
To  consider  the  rotor  effects  it  should  be  remembered  that    u)~ 
J\.fj\  here  K  is  the  relative  motion  of  a  rotor  on  the  body.    Thus 

for  a  rotor  there  will  be  a  combination  of  terms  considering  the  absoiul 
angular  velocity  of  the  OXYZ  system  plus  the  angular  velocity  of  the 
rotor  with  respect  to  the  body. 

Thu .       „ 

s  Ko  +U)  x  Ho     where  U)  »J.L*  r\_ 

dt 

It  should  be  noted,  that  for  a  symmetrical  rotor  or  disk  the 
product  of  inertia  terms  are  zero  where  the  disk  is  at  least  a  three 
bladed  propeller  .    The  two  bladed  propeller  has  a  product  of  inertia 
dependent  on  its  angular  position.    This  case  will  not  be  considered, 
but  the  three  or  more  bladed  propeller  situation  will  be  used  in  all 
considerations.    Hence  with  a  dynamically  and  statically  balanced 
rotor    /t      =  0  and  the  product  of  inertia  is  zero. 

The  moment  of  inertia  of  the  disk  must  be  considered  from  the 

axis  of  the  rigid  body.    Therefore,  the  moment  of  inertia  must  include 

2 

the  md    consideration. 
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AA 


The  moment  of  inertia  terms  for  the  rotors  will  be  thus  define!  in  the 

most  general  sense 

a  =  Ix  rotor  =  Ix  +  md,    and  aligned  with  the  reference  axis 

b  -  Iy  rotor  =  Iy  +  md  2  and  aligned  with  the  reference  axis 

A  A  2 

c  =  Iz  rotor  =  Iz  +  nrcL  and  aligned  with  the  reference  axis 


->       -^       ->  -==> 

G  =   dH  ■  dHo     +  dHrel 

<3FE     dt         dt~"~~ 


Recalling  the  derivation  of  moments  for  the  rigid  body  without  the 
rotor,  the  same  method  applies. 
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fcr  " 


Combining  this  to  a  system  consisting  of  a  rigid  body  plus  a  rotor 
the  following  results  occur 

L 


-» 

* 


c    |>a3  -    3bcZil1  +  IyA3il?  +  Jx/atX2  *  *x*A2fQ 


-*    ' 


*l  *  H?  +telf,  + 


k 


jrti        A?        A3 


Combining  and  using  aerodynamic  notation. 

->  -*r  *     -      t 

s      LAI  '  "7:  *  alfi  ;'B)     ,;   4  c1"^  -  bi^R 

+  T[bQ  -  bKo  -  (A-C)Ht  -  E(P2  -  R?)  +  afc-jR  .  c|^  ^ 

+  k[cR       EP  *  off  ?  +  (B-A)PQ   <   EQE  +  b»f?P  -  9Jf    p] 


This  result  concurs  with  Etkin,  page  116. 

Therefore ,  using  a  standard  aerodynamic  notation  the  moment: 

G=iL+jM+kN 

and  the  scalar  quantities  are  as  follows 
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•    a  K  r  R 


,  '  is,  the  force  and  moment  equations  derived  in  this  appendix 
are  applicable  to  a  rigid  body  of  constant  mass  where  the  body  is  in 
motion.    It  also  considers  the  effects  of  a  three  bladed  or  more  rotor 
system  whicl    Ls    attached  on  some  fixed  point  of  the  rigid  body,,    The 
rotor  effects  may  be  included  or  discarded  depending  on  the  magnit* 
of  rotor  inertia  and  the  relative  velocity  of  the  rotor . 

It  should  be  noted  that  the  derivations  are  done  in  a  general 
form,  but  in  the  concluding  equations,  the  usual  aerodynamic  notation 
is  substituted  to  bring  the  equations  to  this  specific  problem  ,    Further- 
more „  missile  or  aircraft  shape  and  density  has  been  considered  sym- 
metrical about  the  xz  plane »    Thus,  the  product  of  inertia  terms  Ixy  and 
Iyz  have  been  considered  zero.    No  generality  is  lost  by  this  and  if, 
for  some  reason ,  the  body  is  sufficiently  unsymmetrical ,  the  terms  < 

be  quickly  added  in0 

4 
The  following  is  taken  from  Chapter  4  of  Et   in    .    The  concept 

otion  on  the  control  surfaces  is  shown  in  the  following 


<- 


■  Hs 


c 


86 


'1  e  LaGran  ,e  equation  of  notion  in  a  moving  reference  frame  is 


I 


*'        a  "K 


T  =  kinetic  energy  of  the  system  relative  to  the  chosen  frame  of 
reference 

J^  -  -3_1  =  generalized  force 
3  ZK 

W  =  work  done  on  the  system  by  the  external  forces  which  act  upon 
v 

q.h~  generalized  coordinate 

H  =  aerodynamic  moments 

<s  -  rotation  of  control  surface  in  radiar 

Kinetic  energy    T=       I J  5  j_T  j.g 

9  \ 

For  elevator  motion 

Fe  =  generalized  elevator  control  force 
CJ  Wi  =  work  done  by  inertia  forces 
Work  done  by  the  control  forces  is 
cf(<fe)  -  PJ        ♦   Tj 

> 

ds .    "  o  or  -  k  to  control  siirface 

h 

itage  of  hydraulic  cylinder  to  - 
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Evaluating     &  Wi 

consider  dFi  =  dm 

dt 


recalling  from  previous  derivation  that  . 

dt         dt 


Then 
dFxl  = 


dV  8  -       Ry  -  x(Q2  +  R2)   +  y  (PQ  -  R)   +  z  (PR- 

dt 


dFyi  and  ..1  z  are  similar  in  for-  . 

Looking  at  an  elevator  assembly,  it  can  be  seen  that  assumin 

lamina  in  the  xy~plane,  the  displacement  in  the  C^  direction  is 

amount      S  ( i 

The  work   .one  by  the  inertia  forces  is  £-yj    m  \  f*     c  (  S  r 

£ a  fdVoz  /°e  dm  -  (FR-Q)  1x7°  edm  ~  Rr   +  P/V  A?n 

jtj^t  y  dt  x  / 

The  coriolis  terms  2  Py   and  -2  Qx  are  zero  as  there  is  no  work  done 


in  those  directions.    Furthermore  z  =  0  because  of  the  assumed  lamin 

exeac. 


'Yoz  j/°  e  dm  -   m^l^p 


dt  "     z 

rPP   -  Q)  j xf*     'ri   e  Pe      (PR  -  5) 

A 

P^  :  i   cf   th 
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dm    '■   0 


By  symmetry  /  P 

dfe  '^  "  ?    ■  (**-£•)    '  He  +  ; 


F 


j  ,       Pex   (PR-Q)   -  He  4  Fe         (Ktkin  U.8,13) 


Similarly  tor  the  rudder 
I 


tS      -        -\-\     -  (R  +  P°)   P^  -  (RO  -  P)   P      *  H_  +  F 

1  .  >  * 


and  for  the  iileron 
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APPENDIX  II 


ROLL,  YAW  AND  PITCH 


The  derivation  of  the  yaw,  roll  and  pitch  rates  is  given  here  to 
show  the  actual  motion  taking  place.    Frequently  the  motion  is  con- 
sidered small  and  small  angle  approximations  are  used.    These  will  be 
noted  at  the  conclusion. 

First,  consider  the  Eulerian  Axis  used  in  the  aeronautical  field 
shown  below.    The  rotation  takes  place  in  the  order  of  yaw  (    ^  ) , 
pitch  (   O  ) ,  and  roll  (  ft  ) , 


Z.?,  Zt. 


Figure  AII-1      EULERIAN  AXIS  SYSTEM  SHOWING  ROTATION  ABOUT  THE 
Zs  Ys  AND  X  AXIS 
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The  Initial  conditions  at  X    ,  Y    ,  and  Z     and  the  relations  of 

o       o  o 

m  .  -     c  in  be  related  to  P,  Q  and  R  of  the  body  axis  at 
x,  yf  z.    Jn  the    ibove  dia  p  m   consider  the  angular  velocity  P  which 
is  about  the  OX  axis.    It  must  be  the  sum  of  the  velocities  about  the 
OZ    ,  OY.  ,  and  OX.,  axis  respectively. 

OZQ  makes  an  angle  —£    +    3     with  OX.    The  rate  thus 
V  Cos  fflX  h=  A)  -  -  *P  Sin  0°    °Y    is  perpendicular  with  OX  so 


0 


cos  =0  and  OX  is  parallel  with  OX    or  the  motion 

2 


^   y 


Thus  the  summation  P  =  P  -  (p    sin  G 


To  evaluate  Q  the  angle  YOZ     mus  t  be  determined .    Specifically  the 
value  f     cos  YOZ     is  desired. 

cos  y^  os  ybz0  cos  zCcfe. 

O  £       © 


■  cos  (TT  -  0 )  cos  e 

85  sin  fiS  cos  0 


us  the  contribution  of  OZ     axis  motion  is    (f}    sin  p     cos  ©>     . 
The  angle  between  OY  and  OY    is    f>   so  the  motion  here  results  as 

0 

cos  .    OY  is  perpendicular  to  OX„  hence  the  angular  motion 

cos      iL      =  0. 

Z 

•o 

Thus  Q  =     9  cos    j£  +    ^  sin  ^   cos  © 

In  the  third  angular  velocity  motion  R  the  yaw  motion  appears  as 

<&  OB 

cos  ^    cos   &  while  the  motion  of  pitch  is  B  cos  (^  *  0)  —  8  ^ 
The  rolling  term  vanishes    0  cos    -HI       =  0. 
Thus  R  =  -        -    .in  ^    +    *P  cos  ^    cos   O 


Hence  the  relationship  of  the  body  yaw,  pitch  and  roll  rates 
(P,  Q,,    and  R)  to  the  designated  reference  system  angular  velocities 
*P    ,    G  and    0  is: 

F ■  -    <f>  -    qJ  sin   6 
Q  =     0  cos  p  +    $    sin    <£>  cos  © 
R  = ■     lj)  cos  <Z$   cos  O    -  Q   sin  0 
These  results  concur  with  Duncan  page  81  and  Etkin  page  116. 
It  should  be  noted  that  for  small  angles  the  approximation  belovv 
is  often  used: 

q  =    e 

R=      <P 

A  second  consideration  which  must  be  made  regarding  movement 
about  the  reference  axis  is  the  velocity  or  displacement  with  respect  to 
time.    It  is  desired  to  obtain  the  velocity  with  respect  to  the  fixed  fi  a  e 
or  inertial  axis  i*  x,  j*y   ind  - 
Consider  first  that  in  scalar  values 

dx*  =  U 

dT~       1 

dyj!  «V_ 
dt  l 

dz*  =  W, 

It  should  be  noted,  however,  that  in  the  three  position  system  shown 
on  the  previous  page  that 


). 


u    =  u 

3 


V     =    ,  cos  p    -  W     in     <fi 
W  ,  =       sin  ft  +  W  cos  ^ 


in  :         W,    =  -U«  sin  <3    ♦  W-  cos  O 


V3 

U?  =  U.  cos  A   +  WV,  sin  O 


and         U,   =  U«  cos  if 

sin  ^ 

Vx    -  U2  sin   f 

i    y      COS   ^ 

w,  =  w 

Combining  the  1 ,  2  and  3 

position  terms 

dx*  =  U  cos  O  cos  ^  +  V  (sin  p$   sin  G  cos  ^  -  cos  0  sin  ^  ) 
.it 

+  W  (cos  $  sin  ©  cos  If  +  sin  <jf>   sin  <f  ) 

dy^  =  U  cos  &  sin  <f  +  Y  (sin  £$  sin  ^  sin  if  +  cos^S    cos  ^  ) 
dt 

+  W  (cosp*  sine   sin  f  -  sin  <fi  cos  tf>  ) 

dz*  =  -U  sin  e  +  V  sin^f   cos  e  +  W  cos  &>  cos  ^ 

dt 

Thus ,  the  velocity  or  rate  of  change  with  respect  to  the  inertial 
or  earth's  reference  axis  in  this  case  is  expressed  in  these  terms  as 
derived.    These  velocity  terr  s     re  derived  in  Etkin  Page  102  and  coincide 
with  the  derivation  shown. 

Therefore,  this  appendix  accounts  for  the  changes  in  angular  and 
directional  displacement  of  a  rigid  body  with  respect  to  a  fixed  reference 
axis  sy:tem. 


APPENDIX  III 


NOTES  ON  DETERMINANT  MANIPULATION 


The  body  of  this  project  contains  several  determinant  expansions 
where  the  order  is  doubled,     Chese  manipulations  do  not  alter  the  values 
of  the  determinant  and  the  resulting  values  of  the  performance  functions 
are  the  same. 


The  simple  deter;  Ln   n1 


Al 


a; 


A3 


B3 


CI 


( 


C    ■; 


A 


evaluated  will  yield  the  following  result 
2C3  +  B1C2A3  +  CLA" 
C     B2A      •    AIG2B3  -  BIA2C3J 
This  expansion  by  the  equations  -1+1  =   0  will  result  in  the  following 
determinant . 


1 

0 

B] 

0 

c 

0 

1 

1 

0 

0 

0 

0 

A2 

0 

C2 

0 

0 

0 

-1 

] 

0 

0 

0 

A 

■ 

B3 

(.   . 

0 

0 

0 

0 

0 

-.1 

1 

The  resulting  evaluation  will  yield  the  same  value  as  result  (A3-1K 
Furthermore,  if  the  terms  below  the  main  diagonal  are  allowed  to 
remain  in  their  original  colui  n  3  Ln  lieu  of   /hating  one  place  to  the 


\ 


right  s    .one  in  the  nt  dn  body,  the  resulting      /|     is  -itill  the 

ihe  results  <  ■  cofactors  for  particular  performance 

functions  ren  ain  unchanged 0    For  example  it  is  desired  to  evaluate 


cefaclor 


of  the  oi    fina]     etei     Lnant,    Thus: 


"A3" 
A2 


-  -Ci 


1 

BIS 
i 


C2        =*(-!)     /A2C3-A3C2/     A      . 


In  the  expanded  determinant  the  equivalent  cofactor  is 

1 
aa_  -  _o_  -  M  -  —  &■  —  -€h  —  -0" 


■1 


0 


0 


A2        B2 


A3 


C2 


C" 


- 


I 

I 

0 

0 

A2 

G 

A2 

0 

C2 

= 

0 

1 

0 

0 

0 

1 

0 

A3 

0 

0  A3         0  C3| 

.)  JklGZ  -  A3C27 

The  result  of  Equation  A3-3  equals  Equation  A3-2  ,    With  the  v  dues  below 
the  main  diagonal  left  in  their  original  columns,  the  following  occui 


13 


1/% 
*3 


M- 


B 


a; 

o 

0 

0 

a; 

0 

0 

D 

\Q 


Bi2 


-i 


D 


•C4- 

C2 
i. 
C  3 
-] 


-0- 

0 
0 
0 

0 


By  inspection  of  columns  two  and  six,  two  reductions  are  made  to  tru 
following* 

\2  0  C: 

I  C 


13 


A, 


f: 


C: 


=  f-V 


/A2C3  -  A3C2/ 


A        1 


The  resulting  equation  (A3-4)  is  the  same  as  equation  (A3- 2) . 

A  similar  examination  using  a  four  by  four  determinant  system 
was  made.    This  was  done  to  insure  that  the  simplified  effects  oi  a  three 
by  three  or  lower  order  system  did  not  cover  up  an  error  in  a  four  by  four 
or  higher  order  determinant „    The  same  manipulations  can  be  performed 
and  the  resultant  value  of  the  determinant  is  not  altered . 

Thus,  the  individual  has  the  option  of  placing  his  values  in 
either  column,,    The  location  of  the  values  is  dependent  on  desired 
location  of  the  pick-off  point  on  the  block  diagram.    The  final  result 
is  not  altered  c 
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